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ABSTRACT 


MRP system has to be used effectively to reduce the costs and due date 
performance. One of the unaddressed main problem in MRP is the 
uncoordinated arrival of different subassemblies at an assembly centre. 

In this study, this problem has been addressed. A multi-stage, 
multi-product hypothetical production system is simulated to generate data 
about uncoordination. Then, a heuristic which seeks to achieve a better 
coordinated scheduling of all parts is developed. 

Preliminary research indicates that such a heuristic improves the 
performance in a significant manner. Particularly, the aggregate lateness 
of all end items and the actual inventory carrying costs is reduced 


substantially 
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CHAPTER 1. INTRODUCTION 


Owing to the pressure created by the economic environment and loss of 
markets to competition f rom countries throughout the world there is a 

renewed interest in the theory and practice of material management. High 

interest rates and disinflation have intensified the strategic role of 

material management. Additionally, the importance of material management, 
more often referred to as inventory management, for providing better 
customer service in terms of minimizing production lead times, responding to 
the market demand chauiges and meeting customer delivery dates, producing 
products at least cost has been heightened by the embarrassing success of 
f oreign competitors. 

Seventy percent of the costs involved in producing a product is on 
materials. Hence, Inventory management has become the frontier for 

maximizing profits in 1990s [15]. In the area of manufacturing inventory 
management, the most successful innovations are embodied in what has become 
known as material requirement planning, popularly referred as MRP. MRP has 
become a new way of life in production and inventory management, displaying 
older methods in general [16], Inventory management is incomplete unless an 
MRP system is included as em integral part of inventory management. 

MRP is a Production & Operation Management computer information system 
that dynamically develops a schedule of planned orders of materials in each 
time period of the planning horizon [15]. Through a set of logically 
related records, procedures, and rules for decision making, MRP translates a 
master production schedule ( MPS ) into time-phased net requirements for 
each inventory item needed to implement the schedule and spells out the plan 
for covering these requirements. The system specifies materials to be 
procured and the actions to be taken, and determines when they are needed 
during the planning period. 



With MRP, requirements for large raw materials and work-in-process 
inventories are drastically reduced, and better control exists over 
quamtities and timing of deliveries of raw materials, sub assemblies, and 
assemblies to production. These controls help to reduce labour, 
inventories, and overhead costs. 

1.1. Scope of the present work 

MRP has become a popular buzz-word which has impressed business 
associates. MRP is used throughout industry. But, to be in the competitive 
edge over other companies MRP has to be used effectively. It has to be 
modified to address the ground realities of the environment and to keep the 
capital blocked in the shop floor as work-in-process inventory as low as 
possible. 

It is the contention of Dr. Eliyahu Goldratt, f ounder and inventor of 
optimized production technology ( OPT ), that the management of a 
manufacturing organization is made easier if efforts are directed toward the 
control of all critically constraining resource ( CCR ). On shop floor of 
multi-stage, multi-product system using MRP, these CXT^s are coordination of 
all levels of sub assemblies, components, and parts. 

This motivated to study the MRP system and to address this 
uncoordinated scheduling and to find the ways, if any, to reduce 
work-in-process inventory and improve the delivery schedules. The objective 
of this study is as follows; 

1. To develop an MRP/Production system simulation model. 

2. To use this model 

(i) To study imcoordinated behavior of components and parts. 

(ii) To find some ways to attain better coordinated scheduling. 



1.2 Organisation of the Thesis 

Chapter two takes into the world of MRP and gives a helicopter view 
over all concepts, assumptions.mathematics of MRP and also shortcomings of 
MRP. 

Chapter three provides a glance of literature survey, in the area of 
MRP, scheduling and multi-product assembly system that was done bef ore 
starting the work. 

Chapter four describes the new heuristics developed and the idea behind 
it. This chapter also discusses the various performance criteria considered 
to evaluate the heuristics. 

Chapter five provides the description of the system being modeled in 
the present study, with associated assumptions and the performance criteria 
used to evaluate the system. 

Chapter six outlines the general features about the implementation of 
simulation model, with some details about the working of different modules. 

Chapter seven describes the study that was perf ormed and the 
statistical analysis carried out on a set of eight sample problems. The 
detailed result of one sample problem and the summarized results of all 
eight problems are tabulated. 

Chapter eight discusses the conclusion arrived from the study and the 


scope for future work. 



CHAPTER 2. AN OVERVIEW OF MRP SYSTEM 


This chapter provides a general overview of Material Requirement 
Planning ( MRP ) system. This chapter mainly deals with the objectives, 
purpose, assumptions, and prerequisites of MRP, its inputs and outputs, and 
its role in manufacturing, planning & control system. We briefly review its 
working, mathematical models available for analysis of the decisions to be 
taken that greatly affect the system performance. 

MRP is a formal, computerized information system that integrates the 
scheduling and control of materials. Through a set of logically related 
records, procedures, and rules for decision making, MRP translates a master 
production schedule (MPS) into time-phased net requirements for each 
inventory item needed to implement the schedule and spells out the plan for 
covering these requirements. The system specifies materials to be procured 
and the actions to be taken, and determines when they are needed during the 
planning period. [15] 

2.1. MRP objectives 

MRP is necessary because of the volume of materials, supplies, and 
components involved in producing a company’s product line and the speed with 
which management needs to react to the constant dynamic changes of the 
economy. An MRP system has the following objectives: [15] 

* Improve customer service, thereby increasing sales and lowering 
prices 

* Reduce inventories 

* Provide a change-sensitive, reactive manufacturing system 



* Improve efficiency through: 

* Reduced idle time. 

* Reduced setup costs. 

* Avoidance of unplanned delays. 

* Fewer camcelled orders and changes in order quantities. 

* Information for schedule planning before actual release of 
orders. 

• Aiding in capacity planning. 

2.2. Purpose of the MRP system 

The main purpose of an MRP is to control the inventory levels, assign 
appropriate priorities, and plan capacity of the manufacturing system. 
Controlling inventory means ordering the right part at the right quantity 
and time. While assigning priorities among different assemblies it is 
required to consider due dates and inventory carrying costs. Capacity 
planning is to decide on machine capacities with a reasonable capacity 
utilization for a schedule of due dates and throughput. Thus an MRP system 
is designed to answer the question of what need be produced (i.e., what 
capacity is required and what is to be subcontracted) to meet a given master 
production schedule. 

2.3. Assumptions & Prerequisites of MRP 

MRP is effective for companies with the following production 
characteristics of manufacturing operations: 

* Products are produced by using a defined sequence of materials 



(components, parts and sub-assemblies). 

* A dependent demand exists. 

* Demand for components is variable and/or discontinuous in nature. 

Proper use of an MRP is based on following assumptions 
and prerequisites: 

* The final product consists of component parts that are uniquely 
identified. 

* A master production schedule exists, and is stated in 

biil-of -material terms. 

* A computerised bill-of -materials exists at planning time. 

* Inventory records provide data on the status of each item. 

* Lead times are known for the individual items. 

* File data are accurate and up-to-date. 

* All components go in to and out of inventory stock. 

* Discrete dibursement and usages of component material is required. 

* Component parts for the final product are needed at the time of 
assembly order release. 

* Independent processing of manufactured items is done. 

The above requirements are critical for setting up an MRP. 
Implementing an MRP requires that MRP data files, shop floor control, and 
data integrity be well established. In addition, the bill-of -materials 
files must be computerised, which is itself a major task. 

2.4.Role of MRP in Manufacturing Planning and Control Systems 

The field of manuf acturing plaiming and control is concerned with 
management problems and techniques for their solution, as well as with the 
linkages or interactions among particular problem areas. This full system 



includes all data inputs, system modules, and feedback connections. The 
front end section of the MFC system is the set of activities and systems for 
the overall direction settings, which produces master production schedule 
(MPS). The back end or execution section of the system deals with the 
detailed scheduling of the factory and with managing materials coming from 
vendor plants. 

MRP is the central system in the engine portion. It has the primary 
purpose of taking a period-by-period (time-phased) set of master production 
schedule requirements and producing a resultant time-phased set of 
components/raw material requirements. In addition to MPS inputs, MRP has 
two other basic inputs. A bill of material shows, for each part number, 
what other part numbers are required as direct component. The second basic 
input to MRP is inventory status. 

The MRP system serves a central role in material planning and control. 
It acts as a translator of the overall plans f or production into the 
detailed individual steps necessary to accomplish those plans. It provides 
information for developing capacity plans, and it links to the systems that 
actually get the production accomplished. Due to this, the MRP process is 
where many companies have chosen to start the development of their formal 
MFC systems. 

2.5. The mechanics of MRP 

Material Requirement Planning is a technique or a set of systematic 
procedures for managing inventory in a manufacturing operation. There are 
four basic steps to the MRP process: 

1. net 


2. lot 



3. offset 


4. explode. 

^ Netting is subtracting of on-hand and on-order quantities, from gross 
requirements. Lotting is the determination of individual batch or order 
quantities for manufacturing based on the calculated time-phased net 
requirements. There are many lot sizing models or techniques which will be 
discussed later. Some of the techniques involve economic balancing of 
set-up and inventory carrying cost and others are simpler rules such as 
using a fixed nximber of period requirements. Offsetting is the 
determination of the appropriate time for release of planned orders so that 
they will be completed in time to satisfy demand. The planned order release 
data is determined by subtracting the lead time from the date of the 
earliest net requirement which it is intended to satisfy. Explosion is the 
calculation of gross requirements for the next lower level components which 
result from planned orders. Elxplosion may take into consideration shrinkage 
or yield of the components in the manufacture of the higher level item and 
the the quantity of each component required to make the higher level item. 
The "quantity per" factors for the components are multiplied by the planned 
order quantities for the parent item and then, if appropriate, the result is 
divided by the yield factor for the component to give the requirements for 
the components . 

The process of determining net requirements, lot sizing, and 
off -setting of planned orders then performed for the lower level components. 
If the lower level components are also manufactured items, their planned 
orders are exploded into gross requirements for their components and so on 
until planned orders are generated for purchased parts and materials. 

In an on-going situation, each MRP planning cycle reevaluates and 



adjusts the plans generated in the previous cycle. It is this ability to 
reevaluate plans in light of changing conditions and signal appropriate 
corrective actions which makes MRP such a valuable tool in manuf acturing 
control. This adjustment capability is an important part of the MRP 

methodology. Frequent adjustment procedures cause frequent changes in 
plants which causes "nervousness" also. 

2,6. Lot sizing in MRP 

In MRP, lot sizing is activity of translating into an order quantity 
the net unfiled requirements for a given item (while ignoring scheduling 
decisions). Order quantity is a question of minimizing the costs of 
ordering ( or setup ) plus the cost of carrying the item in stock. 
Certainly, it is desirable in MRP to generate the least-cost p\irchasing and 
manufacturing orders. 

There are many approaches to selecting the least-cost order quantity, 
which are discussed extensively in the MRP literature.! Orlicky, 1975; 
Berry, 1972; Wagner-Whitin, 1958; Dematteis, 1968; Gorham, 1968; 
Silver-Meal, 1973 ). Among these approaches, Wagner-Whitin algorithm is an 
optimizing technique based on dynamic programming model. It is popularly 
used as it has been proved to be computationally efficient. It is described 
below in brief for the sake of completeness. 

2.6.1. Wagner-Whitin Algorithm 

Wagner-Whitin algorithm 1251 evaluates all possible ways of ordering to 
cover net requirements in each period of the planned horizon. Its objective 
is to arrive at the optimum ordering strategy for the entire net 
requirements schedule. The Wagner-Whitin algorithm is "elegant" in that it 



reaches this objective without actually having to consider, specifically, 
each of the strategies that are possible. The algorithm is as follows : 

Let quantity needed, 

i^ = interest charge per unit of inventory carried forward to 


period t +1; 

s^= ordering ( or setup cost) 
x^= quantity ordered (or manufactured) 
where t denotes period, t = 1,2,...,N. 

It is assumed that all period demands and costs are non-negatives. The 
problem is to find a program >= 0,t = 1,2,..., N, such that all demands 
are met at a minimum total cost; 


Let F(t) = min 


min 


t-i 




t 

I 


h=:J k=h+l 


‘h “k * ‘ > 


+ F(t-l) 


where F(l) = s^and F(0) = 0. That is, the minimum cost for the first t 
periods comprises a setup cost in period j, plus charges for filling demand 
dj^ , k = J + l,...,t, by cauTying inventory from period j.plus the cost of 

adopting an optimal policy in periods 1 through J -1 taken by themselves. 

« » 

The Algorithm at period t , t = 1,2,...,N, may be stated as 

«« mm « 

1. Consider the policies of ordering at period t , t - 1,2,..., t , 

m mm «« m 

and filling demands d^, t = t ,t +l,...,f , by this order. 

m 

2. Determine the total cost of these t different policies by adding 

the ordering and holding costs associated with placing an order at period t, 

mm 

and the cost of acting optimally for periods 1 through t -1 considered by 

themselves. The latter cost has been determined previously in the 

* » 

computation for periods t = 1,2 t -1. 

m 

3. From these t alternatives, select the minimum cost policy for 

m 

periods 1 through t considered independently. 

* m 

4. Proceed to period t + 1 { or stop if t = N ) 



2.7. Sequencing in MRP 


MRP provides component due dates by the expected lead time from each 
parent item. Once the due dates have been determined, it is left to the 
shop floor control (SFC) system to meet these deadlines. The sequencing 

rule is used to select the next job to be processed from a set of jobs 
awaiting service. Sequencing rules are normally used to minimize total 
inventory cost. The nature of the dispatching rule employed influences 
delay costs, in-process and final inventory costs, and setup costs. Due 
date based sequencing rules are the best performers in the multi stage job 
shop [5]. Critical ratio rule (CRR), a member of the family of due date 

based rules, performs well in MRP environment [8]. 

2.7.1. Critical Ratio Rule ( CRR) 

A good discussion of several critical ratio rules ( including some in 
actual use at Black and Decker, Hughes Archoff, and Western electric ) is 

provided by Putnam et al [18]. In its most general form, it is computed as 

[time remaining - work remaining] 

Critical Ratio = 

[work remzdning] 

If CRR value of a job is less than one, then, it is behind the schedule 
and it has to be given high priority. If CRR value is greater than one, 
then, it is ahead of schedule and it can be delayed. Simply, it means, the 
least CRR value, higher the priority and vice versa. 



2.8. Shortcominfs of MRP 


Many of the assumptions made in MRP are unrealistic. In this section, 
we will discuss some of these assumptions, the problems that arise as a 

result of them, and the means for dealing with these problems. 

2.8.1. Uncertainty 

Underlying MRP is the assumption that all required information is 
knownwith certainty. However, uncertainties do exist. The two key sources 
of uncertainty are the forecasts for future sales of the end item and the 
estimation of the production lead times from one level to another. Safety 
stock can be included to protect against the uncertainty of demand. The 
manner in which the uncertainty transmits itself through a complex multi 

level production system is not well understood. For that reason, it is not 

recommended to include independent saf ety stock at all levels of the 
systems. Rather suitable safety levels can be built into the forecasts for 
the end item. These will automatically be transmitted down through the 
system to the lower levels through the explosion calculus. 

2.8.2. Capacity Planning 

Another important issue that is not treated explicitly by MRP is the 
capacity of the production facility. Capacity requirement planning (CRP) is 
the process by which the capacity requirements placed on a work center or 
group of work centers is computed by using the output of the MRP planned 
order releases. If the planned order releases result in an infeasible 

requirements schedule, there are several possible actions. One is to 



schedule overtime at the bottleneck locations. Another is to revise the MPS 


so that the planned order releases at lower levels can be achieved with the 
current system capacity. This is clearly a cumbersome way to solve the 
problem, requiring an iterative trial and error process between the CJIP and 
the MRP. 

2.8.3. Rolling Horizon and System Nervousness 

In practice, production planning environment is djmamic. The MRP 
system may have to be rerun each period and the production decisions are 
reevaluated. Often it is the case that only the lot sizing decisions for 
the current plaiming period need to be implemented. Rolling horizon is used 
to refer the situation in which only the first-period decision of an 
N-period problem is implemented. The full N-period problem is rerun each 
period to determine a few first-period decision. 

Another common problem that results when using MRP is nervousness. The 
term was coined by Steele [23], who used it to refer to the changes that can 
occur in a schedule when the horizon is moved forward one period. Some of 
the causes of nervousness include imanticipated chamges in the MPS because 
of updated forecasts, late deliveries of raw materials, failure of key 
equipment or absenteeism of key personnel, and xmpredictable yields. 

2.9. Factors affecting the computation of requirements 

The computation of requirements is complicated by six factors: 

1. The structure of the product, containing several manufacturing 
levels of materials, component parts, and subassemblies 

2. LOT sizing, i.e.,the ordering of inventory items in quantities 
exceeding net requirements, for reasons of economy or convenience 



3. The different individual lead times of inventory items that make up 
the product 

4. The timing of the end-item requirements across a planning horizon 
of, typically, a year’s span or longer, and the recurrence of these 
requirements within such a time span 

5. Multiple requirements for an inventory item due to its so-called 
commonalty, i.e., usage in the manufacture of a number of other items 

6. Multiple requirements for an inventory item due to its recurrence on 
several levels of a given end item. 



CHAPTER 3. PREVIOUS RESEARCH 


The recent growth in the use of material requirement planning (MRP) 
systems has resulted in increased interest in the topic of decision making 
in multi-stage, multi-product systems. MRP has been subjected to extensive 
research. Particularly, there is hardly a scarcity of research in the area 
of lot sizing and dispatching rules, their interaction and their impact of 
MRP system. 

Biggs.J.R.[4] reported that the effect of lot sizing and sequencing 
rules and the interaction between these is significant at the level of 0.01 
significance. He has argued that a particular sequencing rule may work in 
opposite direction to particular lot sizing rule or vice versa. 

Blackburn. J.D., and Millan.R.A. [7] have modeled a multi-stage problem 
analytically to indicate the potential errors inherent in the commonly 
purposed single-pass, stage-by-stage approaches. Then, based on this 
analysis, several simple cost modifications are suggested. 

Biggs.J.R.[5] has studied the effect of using various priority 
scheduling rules to capacity management and shop floor control in a 
multi-stage, multi-product system which uses MRP. He foimd that Critical 
Ratio Rule (CRR) performs consistently well over different levels of 
capacity. 

Grasso and Taylor [11] studied main effects and interaction effects 
among four factors, namely, lead time, buffering alternatives, lot size rule 
and cost value on total cost and reported that maun effects of each factor 
and the interaction between lead time distribution and lot size auid that 
between buffering alternative and cost value on total cost are significant. 

South and Steward [22] report their experience at Hix corporation 
wherein fifty percent increase in lot sizes and planned manufactxired lead 
time resulted in low system tardiness without increase in WIP inventory 



costs. 

As already discussed. Lot sizing, sequencing, determination of planned 
lead times and capacity planning are important decisions in the eff ective 
use of MRP. There is a lot of literature available on scheduling. 

Panwalkar.S.S., Iswander.W. [17] have classified over 100 scheduling 
rules and made an attempt to explain the general idea behind different 
rules. 

Blackstone .Philips and Hoggs [8] have discussed the state of the art 
in the study of dispatching rules, listed 34 dispatching rules and compared 
them. They have concluded that no single dispatching rule developed so far 
consistently produce lower total cost than all other rules under a variety 
of shop configurations and operating conditions. 

Goodwin, J.S., and Goodwin, J.C., [13] have shown that the performance of 
various sequencing rules do not generalize the multi-stage job shop. 

In a recent study, Russell and Taylor [20] have evaluated the 
performance of various dispatching rules in two different product 
structures. The authors constructed a BOM representing a tall product 
structure having four levels and a BOM representing a flat product structure 
having three levels. They concluded that as the BOM gets taller, the 
probability that the job will finish after its due date increases. 

In another study Goodwin and Weeks [13] have studied the effect of 
sequencing rule and frequency of priority updating. They concluded that due 
date based sequencing rules were the best performers in the multi-stage job 
shop. 

Fry, Philipoom and Markland [12] have evaluatedthe performance of 
selected dispatching rules in an imbalanced multi-stage job shop to 
determine whether the performance is consistent for jobs which are rated 
through the bottleneck and f or jobs which by-pass the bottleneck. They 
concluded that sequencing rules which are based on job due date tends to be 



the overall better performers. They also concluded that the results for 
jobs which pass through the bottleneck were not the same for jobs which 
by-pass the bottleneck. 

Potty. V.S. [18] has reported that the choice of a lot size rule has 
significant impact on the total costs and the choice of sequencing rules has 
significant effect on number of stock outs and unit stock outs. He also 
reported that lot sizing using Wagner-Whitin algorithm and sequencing using 
CRR rule does well in MRP environment. 

In a recent study, Sharma.R.R.K., and Potty. V.S. I211 have presented a 
multi-pass heuristic which improves the performance on the number of stock 
outs by an iterative interaction with the lot size rule under fixed 
capacity. 

But, Sequencing in a multi-stage, multi-product system requires a 
coordinated flow of materials throughout the various stages to assure 
components needed for assembly are available when planned. Coordination is 
vital since shortage of a singie item can halt the assembly process, 
resulting in sharply diminishing productivity and increased inventory costs. 

But, surprisingly, no work has been done in this area to question the 


uncoordinated scheduling. 



CHAPTER 4 . HEURISTICS DEVELOPED 


Normally, single pass heuristics are used in MRP systems for 
simplicity. That is, a single lot sizing procedure is applied to determine 
the lots and next a single scheduling heuristic is applied. This procedure 
is popular due to its simplicity. But such a procedure lacks communication 
between the lot sizing module and scheduling module in various respects. 
Biggsl4] has established that in fact lot sizing and sequencing rules 
interact. Hence, there may be scope to improve the solutions by using multi 
pass hexiristics. One such effort can be foimd in the work by Sharma and 
Potty[21], where the inventory holding costs were increased for items or 
assemblies which spends substantially more time in the shop than the 
theoretical holding costs. 

Here, we note that . possibility exists to improve lateness/tardiness 
performance if variance of the arrival times of different subassemblies at 
an assembly centre can be reduced. This means, in other words, that efforts 
made to improve the co-ordinated arrival of jobs at various assembly centres 
could be fruitful. With this preliminary hypothesis we developed a 
heuristic which seeks to achieve a better co-ordinated arrival at assembly 
centres. This heuristic is described below. 

4.1. Heuristic 1 

This heuristic is based on the following idea. 

i) . Find the components which arrived very late and forced 
assemblies to wait. Ebcpedite it by increasing its criticallity. 

ii) . Find the components which arrived very early and waited for a 



long time. Delay it by decreasing the criticallity. 

The steps of heuristics are described below. 

Let T = length of the planning horizon. 

BestSol = Best solution arrived at w.r.t aggregate lateness. 

R = set of components which has to be expedited. 

D = set of components which has to be delayed. 

1. Initialization. Set E = {};D = {}; 

BestSol = a ; 

2. Read input data. 

3. Simulate MRP system for the time horizon considered. 

4. Find the components which have highest aggregate waiting time and 
let j be that component. E = E + {j} ; 

5. Find the components which have less aggregate wadting time auid let k 
be that component. D = D + {k} ; 

6. If the current solution is better than the BestSol w.r.t aggregate 
lateness, then update BestSol = current solution. 

7. Increase the criticallity of all components in set E in steps of 

0 . 1 . 

8. Decrease the criticallity of all components in set D in steps of 

0 . 1 . 

9. If the n^mlber of iterations considered is less than twenty then go 
to step 3 else print final results. 

An increase in the criticallity of a component by increasing the 
CRR-value, it gets high priority and arrives at early period. The opposite 
is the case for decreasing CRR-value. So, the aggregate lateness is reduced 
Md the coordination is expected. Actual inventory carrying cost is also 
expected to reduce. For the sake of comparing the performance of the above 
heuristics with other multi pass heuristics, we reproduce a heuristic due to 



Sharma and Potty with slight modification. 


4.2. Hexiristic 2 

Let R be the set of parts with difference in actual and theoretical 
holding costs exceeding 10 7, 

Let S be the set of parts already considered twice in the iteration 
procedure. 

Let K( j ) denotes the number of times the part "j" has been considered 
in the iterative procedure. 

Let BestSol denote the best solution arrived w.r.t the agg regate 
lateness. 

The steps of the heuristic are described below. 

1. Initialize the following 

S = { } ; K{j) = 0 for all j . 

BestSol = a ; 

2. Read input data 

3. Simulate the MRP system for the time horizon considered to prepare 
the set R. 

4. Select part "j" from set R which has the highest percentage 
difference in theoretical and actual holding costs. Increase the holding 
cost of parts "j" in steps of 0.2 xmtil the lot size changes. 

5. Set K{j) = K(j) + 1; 

6. If ( K(j) >= 2 ), then update S * S + < j } ; 

7. If the current solution is better than the BestSol w.r.t aggregate 
lateness than BestSol 4= Current solution. 

8. R <= R - S . 

9. If R = { } then stop else go to step 3. 



An increase in the holding costs is likely to reduce the lot size and 
reduce the time taken to process that lot. That means, it arrives to the 
assembly shop eau-ly and get assembled. It is expected that it may reduce 
the aggregate lateness of all assembled items and the actual inventory 
costs. But, in the same time , an excessive decrease in lot sizes may 
increase the number of setups, and adversely affect the no of stock outs and 
the total costs. '' 

•4.3. Heuristic 3 

In heuristic 2, there is no effort made to secure co-ordination. It is 
also possible to reduce the aggregate lateness by 

i) . Changing the lot size and arrives the best solution using the 
Heuristic 2. 

ii) . Then, Finding the components which created the uncoordination 
and applying Heuristic 1. 

It is expected that simultaneous changes in the lot sizes and the 
priorities may reduce the aggregate lateness and the components are expected 
to arrive with coordination with one another. 

The steps of heuristicis described below. 

1. Initialize in the same way as in the case of Heuristic 1 and 
Heuristic 2. 

2. Simulate MRP system 

3. Apply the Heuristic 2 and get BestSol. 

4. Starting with this IIBestSol , Apply the heuristic 1 and get updated 


Solution. 
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4,4, Performance criteria 

The performance of the system is analysed using the following 
representative criteria : 

1. Aggregate lateness of all end items, 

2. Total number of stock outs, final products, 

3. Inventory costs, total system. 

4. Totad n\imber of setups, total system. 

Aggregate lateness of end item is considered as a performance criteria 
because this lateness is due to the uncoordinated arrival of all the 
components. 

When the manufacturing organization are facing a tough competition, 
their perf ormance is judged on the basis of the punctuality of delivering 
the pairts as ordered. Number of stock outs is considered as a performance 
criteria. 

Inventory carrying costs indicate the extent to which capital has been 
locked up and has a direct influence on the profits earned. 

The number of setups is also important, especially when the change over 
times from one setup to another setup are significant. 

The Simulation model used to evaluate the above heuristics is followed 


in the next chapter. 



CHAPTER 5. SIMULATION MODEL 


A computer simulation model of a hypothetical multi-stage, 
multi-product production system is used to generate data about the 
uncoordinated scheduling. This chapter provides the description of the 
system being modelled in the present study, with associated 
assumptions. Simulation is defined as " the process of creating the essence 
of reality without ever actually attaining the reality itself". Simulation 
is selected because 

1. The system behaviour under disturbance can be predicted easily. 

2. Detailed analysis of a complex material requirement planning system 
is relatively easy. 

3. Simulation provides flexibility in experimental design. 

5.1. Logic of the model 

The simulation model implements standemd MRP logicC Fig 5.1 ). 

Although a simulation model, it operates essentially deterministically. 
Since product structures and routing are fixed, it processes entities in a 
known manner. Demand by time period for end items are assumed known ( but 
not constant) over a finite horiion. 

5.2. Assumptions 

The system has been designed under the following assumptions: 

1. Demand is deterministic. 

2. Man pmwer is of uniform ability. 

3. Breakdown times and repair times of the machines are deterministic. 














4. Pre-empting is not allowed. 

5. Machine setup and operation times are known. 

6. No work center can perform more thain one operation at any time. 

7. Each lot size, once started must be completed if the component parts 
are available. 

5.3. System description 

The entire system is divided into the following five basic modules: 

1. Input data module. 

2. Planning module. 

3. Shop floor control module. 

4. Simulation module. 

5. Output module. 

Input data module: 

In this phase, all the input data are read, like, the product structure 
from Bill-of -material files, external forecasted demand for the end items 
from the master production schedule file and the inventory status of all 
items from inventory status file. An inventory status file contains data 
for each item in inventory, like, scheduled receipts, ordering cost, 
inventory carrying cost over planning horizon, lead time and quantity on 
hand 

Planning module : 

In this module, a detailed exploration of the end items to all 
subassemblies, components and parts takes place. Then, netting, lot sizing 
and offsetting takes place. The output of this module is the release of 
planned orders to the shop floor to start the production operation and to 
the pxirchase department to plan purchase orders. 



Shop floor control module: 

In this module, all the planned orders for a particular period are 
undertaken for the production. Each part, then, follows their routing 
pattern over various processing centers. Once the operation gets over in a 
particular machine, it goes to next machine and waits in the queue. Once 
the machine is free, the next critical job in the queue will be tadcen for 
the processing. Sequencing rule is used to find the criticality of a Job. 
After all processing is over, the transformed part follows to assembly shop 
, if it has any parent. . Once all necessary components have arrived, the 
assembly operation starts. 

Simulation module: 

This module perf orms the simulation work. It runs and controls the 
shop floorover the planning horizon considered. The planning horizon 
consist of forty time buckets. 

Output module: 

Once the assembly of end item was completed, the end item will be 
dispatched to the customer. This module provides all the statistics like 
lateness, number of stock outs, setup costs, inventory carrying cost and 
total cost incurred over the planning horizon. 

5.4. Size of the shop 

When Baker and Dzielinski [2] tested some of the dispatching rules in 
shops of various sizes, they foimd that size of the shop does not affect the 
relative performance of the rule. This finding is convenient in that it is 
cheaper to simulate a small shop. Buffa [91 concludes that: " Since shop 
size has never appeared as a major variable, it seems that we may be able to 
experiment with small shops and generalize the resulting conclusion". The 



production system used in this study consists of fifteen processing centers 
and five assembly centers. But, it can be varied. 

5.5. Ruies considered 

Three heuristics discussed in earlier chapter are implemented in the 
simulation model. Lot sizing considered is Wagner-Whitin algorithm and 
sequencing rule considered is critical ratio rule. 

5.6. Modeling approach 

Discrete event simulation with next event scheduling approach is used 
in this model. In this approach there is no change in the system status 
between any two consecutive events. The clock time of simulation and the 
system status are updated at the occurrence of each of the events occurring 
in the system. In the model that has been developed, the following are the 
events that causes a change in the system status: 

1. Order release at the starting of the period. 

2. Processing machine loaded/imloaded. 

3. Assembly machine loaded/unloaded. 



CHAPTER 6. IMPLEMENTATION OF THE MODEL 


The heuristics developed in the previous chapter have been implemented 
in an MRP environment simulated on a personal computer. This chapter 
outlines the general features about the implementation, with some details 
about the working of different modules. 

6.1. A word about the implementation 

The essential components of the modelsuch as a clock mechanism, data 
structures to represent transactions, resources, and queues and collection 
of essential statistics were identified. job. To cope up with this lengthy 
process of building a simulation model, a number of special purpose 
simulation languages like Simula, GPSS, Simscript, etc, have been developed 
and widely used. Although all these special purpose languages have passed 
through several versions, none of them cam be conveniently used to implement 
the user defined heuristics. Pascal is a highly structured lamguage and 
often highly sophisticated data structure, and also support dynamic 
allocation of memory. Because of these reasons data structures of linked 
lists type can be created to handle the various events, and thus the memory 
requirements does not depend on the length of the simulation. Because of 
these reasons, Pascal is used for the implementation of the present 
situation model. 

6.2. Working of different module 

In this section, the working of different procedures in our simulation 
model is described. 

6.2.1. Input data module 

In this modules three procedures namely, read_assly_data, 



























read_deinand_data and read_inventory_data are called. read_assly_data reads 
all the details from bill-of -material file. read_deinand_data reads all the 
external demands from demand file for all the assembled items. 
read_inventory_data opens the inventory-status file amd reads the inventory 
status, like scheduled receipts, ordering costs and inventory carrying cost 
over planning horizon , lead time and quantity on hand. 

6.2.2. Planning module 

In this module, there are three procedures, namely, do_planning, 
do_lot_sizing and explode. do_planning cadis the explode. . explode explores 
all the product structures starting from the end items to all purchased 

parts to find gross requirements of all sub-assemblies, parts and 
components. do_planning calls do_lot_sizing after its complete netting. 

do_lot_sizing is the implementation of Wagner-Whitin lot sizing algorithm. 
Then do-planning does the offsetting using the lead times of all components. 

The final output of this module is the detailed planned orders for all 

components over the time horizon. 

6.2.3. Shop floor control module 

This module is the imp>ortant module in the model. It controls all the 
activities of shop floor from receiving the orders from planning module to 
sending the final finished products of all end items to finished goods 
inventory store. 

The working of different procedures are described below: 

initialise_the_item accepts the planned orders and does the following: 

i) . If the item needs processing, then, depending on the first 

machine where the first operation has to be done, the item is inserted in 

the queue before that machine using insert. 

ii) . If the item doesn’t need processing, then, the item goes to 

assembly machine on which its parent item is going to get assembled. This 



Job has been done by procedure «earch_for_parent. 

insert is a procedure used to insert a item in the queue bef ore the 
processing machine. In case, the machine is free, the item is loaded 
immediately using load_the_processing_inachine. 

insert_in_assly_mc receives all the orders for the assembled items from 
the planning module and the assembling machines in the assembly shop has 
been informed and asked to keep track of incoming sub-assemblies, components 
and parts. 

search_for_parent is another important procedure in this model. The 
input to the model is either a finished component or parts in the plant 
itself or a purchased component from the outside which does not need any 
processing. The incoming item finds its parent item in the following way. 

i). If no_of_parent is one, then, the item has to be placed before 
the assembly machine where its parent is getting assembled. 

ii). If the item has more than one parent item,then the parent item 
which gets assembled at an early stage gets priority. In case if there are 
more than one parent item getting assembled at the same period, then the one 
which needs less quantity gets priority. 

load the processing machine is a procedure which will be called 
whenever the processing operation on a particular machine getscompleted. 
Then, the finished component is unloaded and updated using update_the_item. 
The next highly critical item is loaded on this machine. The finishing time 
of this machine is updated. 

update_the_item does the following; 

Once a particular job is gets completed, then the job has to be sent to 
the next machine where the next operation has to be carried out. insert is 
called for the above purpose. In case if the item has completed all 
operations, then, the finished item willbe sent to assembly shop using 



search_the_pareiit. 

load_the_assly_mc loads the assembly machine when all the required 
components are available in required quantity. It also updates the 
information like finished time, no of stock outs in case if it is late etc. 

6.2.4. Simulation module 

This is the module which does all the simulation work. 

For a particular period, sequencing initialisation does all the 

necessary initialisations that has to be done at every time 

period. start_the_process starts loading all the processing machine if they 
free. find_next_executable_cvent finds the next executable event that may 
occur in the shop and change the system. update_clock_time updates the 
current clock time to time at which the next executable event occurs. 

do_sequencing runs the model f or the time horizon considered . In this 

study, the time bucket consists of fourty time buckets. 

6.2.5. Output_module 

This module collects all the related statistics using 
collect_and_jprint_statistics. This also prints the complete result 



CHAPTER 7. TEST RUNS fc RESULTS 


This chapter gives the computational results generated and the 
statistical analysis carried out on a set of eight sample problems. 

7.1. Simulation Run 

The simulation was performed over a planning horizon. The planning 
horizon considered was forty periods. The time period considered was a 
week. The statistics were collected only after a preliminary r\in for five 
periods. The sample problems were selected in such a way that it has all 
different product structures like flat, tall and mixed. Figures 7.1.1 to 
7.1.23 are the different product structures considered in oior study. Every 
problem contains three to five different product structures. Table 7.1.1 
gives the details about the product structures contained in eight different 
problems. Inventory carrying costs from one period to another period and 
setup ( ordering ) cost over all periods were considered as variable. 

For illustrative purpose, the details of problem no 4 is given fully. 
The number of end products in this problem is four and total number of parts 
in this problem is forty five. The subassembly details are given in the 
table 7.1.2. The work center where the assembly process takes place and the 
time taken to assemble per unit are given in the table. 

The routing details of the parts are given in the table 7.1.3. The 
numbers within the brackets indicate work centre where the process takes 
place and the time taken to process per unit. Hie machining sequence of a 
part is the same as given in the table. 


7.2. Discussion of Results 
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TABLE 7.1.2 
SUB-ASSEMBLY DETAILS 


PART 

NUMBER 

WORK 

CENTER 

TIME PER 
ASSEMBLY 

1 

16 

40 

2 

18 

45 

3 

19 

30 

4 

17 

25 

5 

20 

15 

1 9 

17 

20 

20 

16 

30 

27 

18 

34 

34 

26 

20 

39 

19 

45 

42 

16 

15 
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Table 7.2.1 
SUMMARY OF RESULTS 
PROBLEM NO 1 



SOLUTION 


ONE SOLN 

PERCENTAGE 

IMPROVEMENT 

HEURISTIC 
TWO SOLN 

PERCENTAGE 
I MPROVEMENT 


THREE SOLN 


I MPROVEMENT 


AGGREGATE 

LATENESS 

ACTUAL 
INVENTORY 
COST 
( ’ 000 ) 

8732.61 

12.07 

8005.35 

7.62 

8.32 

39.54 

8228.02 

9.61 

5.77 

23.39 

j 

8005.0 

5. 25 

8.32 

58 . 3 


TOTAL 

COST 

NO OF 

SET UPS 

( ’ 000) 


42 . 42 

698 

37 . 38 

698 

11.88 

0 

39 . 68 

699 

6 . 48 

-0. 14 

35.31 

619 

16 . 76 

0. 14 








51 


TABLE 7.2.2 
SUMMARY OF RESULTS 
PROBLEM NO 2 



AGGREGATE 

LATENESS 

ACTUAL 
INVENTORY 
COST 
( ’000) 

TOTAL 

COST 

( ’ 000) 

NO OF 
SET UPS 

NO OF 

STOCK OUTS 

STARTING 

SOLUTION 

605.86 

11.82 

28.15 

791 

39 

HEURISTIC 
ONE SOLN 

605.86 

5.81 

22. 13 

791 

39 

PERCENTAGE 

IMPROVEMENT 

0.00 

54.00 

21.38 

0.00 

0.0 

HEURISTIC 
TWO SOLN 

605.56 

5.98 

22.66 

791 

39 

PERCENTAGE 

IMPROVEMENT 

0.0 

50.63 

19.71 

0.0 

0.0 

HEURISTIC 
THREE SOLN 

605.86 

4.54 

20.82 

794 

39 

PERCENTAGE 

IMPROVEMENT 

0.0 

61.5 

25.80 

-0.37 

0.0 


J J 308 ^ 
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TABLE 7.2.3 
SUMMARY OF RESULTS 
PROBLEM NO 3 


STARTING 

SOLUTION 


HEURISTIC 
ONE SOLN 

PERCENTAGE 

IMPROVEMENT 


HEURISTIC 
TWO SOLN 

PERCENTAGE 

IMPROVEMENT 


HEURISTIC 
THREE SOLN 


PERCENTAGE 

IMPROVEMENT 


AGGREGATE 

LATENESS 


9793.73 

8055.33 

17.75 

9793.73 

0.0 

8055.33 

17.75 


ACTUAL 
INVENTORY I 
COST 
( ’ 000 ) 


TOTAL 

COST 

( ’ 000 ) 


NO OF 
SET UPS 


21.88 58.44 801.0 


17.13 44.38 801.0 


21.72 24.0 0.0 


21 . 88 

0.0 


52.47 

10.21 


807.0 

-0.74 


17.13 49.91 801 


21.72 14.59 0.0 


NO OF 
STOCK OUTS 


31 .0 

29.0 

6.45 

31.0 

0.0 

29 

6.45 
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TABLE 7.2.4 
SUMMARY OF RESULTS 
PROBUIM NO 4 . 



STARTING 

SOLUTION 


HEURISTIC 
ONE SOLN 

PERCENTAGE 

IMPROVEMENT 


9. 11 


HEURISTIC 
TWO SOLN 

PERCENTAGE 1.12 
IMPROVEMENT 


HEURISTIC 19976.7 
THREE SOLN 


PERCENTAGE 9.11 
IMPROVEMENT 


ACTUAL 
INVENTORY 
COST 
( ’ 000 ) 

TOTAL 

COST 

{ ’000) 

NO OF 
SET UPS 

92.88 

220.67 

1085.0 

52.01 

144.52 

1085.0 

56.93 

22.04 

0.0 

49.83 

145.95 

1125.0 

46.73 

21 . 29 

-3.6 

43.6 

143.06 

1142.0 

46.57 

22.85 

5.25 
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TABLE 7.2.5 
SUMMARY OF RESULTS 
PROBLEM NO 5 



STARTING 

SOLUTION 


ONE SOLN 

PERCENTAGE 

IMPROVEMENT 

HEURISTIC 
TWO SOLN 

PERCENTAGE 

IMPROVEMENT 


THREE SOLN 


IMPROVEMENT 


AGGREGATE 

LATENESS 

ACTUAL 
INVENTORY 
COST 
( ’000) 

11104.33 

85.95 

10292.84 

47.58 

7.37 

55.43 

11098.40 

49.32 

0.2 

42.62 

10292.84 

32.75 

7.71 

64. 13 


( ’ 000 ) 


NO OF 

NO OF 

SET UPS 

STOCK OUTS 


223.02 758 


20.58 0.0 


228. 12 766 


18.77 I 0.0 
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TABLE 7.2.6 
SUMMARY OF RESULTS 
PROBLEM NO 6 


AGGREGATE] 

LATENESS 


ACTUAL 
INVENTORY 
COST 
( ’ 000 ) 


STARTING 

SOLUTION 


17529.48 33.74 


HEURISTIC 
ONE SOLN 


16980.52 18.45 


PERCENTAGE 

IMPROVEMENT 


3.13 45.3 


HEURISTIC 
TWO SOLN 


PERCENTAGE 

IMPROVEMENT 


17529.48] 

0.0 


33.74 

0.0 


HEURISTIC 
THREE SOLN 


16980.52 18.45 


PERCENTAGE 

IMPROVEMENT 


3. 13 


45.3 


TOTAL 

COST 

NO OF 

SET UPS 

{ ’000) 


76.50 

1457 

61.21 

1457 

19.98 

0.0 

76.69 

1462 

-0.26 

-0.34 . 

61.21 

1457 

19.98 

0 


NO OF 
STOCK OUTS 


66 

68 

3.03 

69 

4.55 

70 

6.06 
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TABLE 7.2.7 
SUMMARY OF RESULTS 
PROBLEM NO 7 



STARTING 

SOLUTION 


ONE SOLN 

PERCENTAGE 

IMPROVEMENT 

HEURISTIC 
TWO SOLN 

PERCENTAGE 

IMPROVEMENT 


THREE SOLN 


AGGREGATE 

LATENESS 

ACTUAL 
INVENTORY 
COST 
( ’000) 

6293.62 

72.48 

6185 

30.72 

1.72 

57.62 

6293.62 

7.86 

0.0 

73.46 

6185 

26.68 

1.72 

63.19 



87.39 


32.33 


66.86 724 


84.02 


65.04 -0.28 
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TABLE 7.2.8 
SUMMARY OF RESULTS 
PROBLEM NO 8 



AGGREGATE 

LATENESS 

ACTUAL 
INVENTORY 
COST 
( ’000) 

TOTAL 

COST 

( ’000) 

NO OF 
SET UPS 

NO OF 

STOCK OUTS 

STARTING 

SOLUTION 

4756.10 

21.83 

89.35 

482 

18 

HEURISTIC 
ONE SOLN 

4493. 10 

8.57 

76.09 

482 

19 

PERCENTAGE 

IMPROVEMENT 

5.53 

60.74 

14.84 

0.0 

5.35 

HEURISTIC 
TWO SOLN 

4756.10 

13.46 

80.94 

498 

18 

PERCENTAGE 

IMPROVEMENT 

0.0 

38.34 

9.41 

3.33 

0.0 

HEURISTIC 
THREE SOLN 

4493. 10 

5.22 

70.77 

4.83 

18 

PERCENTAGE 

IMPROVEMENT 

5.53 

76.12 

20.79 

0.27 

0.0 









PROBLEM 

NO. 

7. IMPROVE- 
MENT OF 
HI 

7. IMPROVE- 
MENT OF 

H2 

7. IMPROVE- 
MENT OF 

H3 

1 

8.32 

5.77 

8 . 32 

2 

0.0 

0.0 

0.0 

3 

17.75 

0.0 

17.75 

4 

9.11 

1 . 1 

9.11 

5 

7.31 

0.2 

0 . 00 

■ 6 

3. 13 

0 . 0 

3 . 13 

7 

1 . 72 

0 . 0 

1 . 72 

8 

5.53 

0.0 

5 . 53 

i 

MEAN 

6.61 

0.88 

5 . 70 

S.D. 

5 . 54 

2.01 

5 . 99 








TABLE 7.2.10 


IMPROVEMENT IN ACTUAL INVENTORY CARRYING COST 


PROBLEM 

NO. 

7. IMPROVE- 
MENT OF 

HI 

7. IMPROVE- 
MENT OF 

H2 

7. IMPROVE- 
MENT OF 

H3 

1 

39.54 

23.37 

58.30 

2 

54.00 

50.63 

61.52 

3 

24.00 

0.00 

21.72 

4 

56.93 

46.73 

46.57 

5 

55.43 

42.62 

64.13 

6 

45.30 

.0 

45 . 30 

7 

52.62 

73.46 

63 . 19 

8 

60.44 

38.34 

76.08 

MEAN 

48.53 

34.39 

54 . 60 

S.D. 

11.95 

25.39 

16.56 











TABLE 7.2.12 


IMPROVEMENT IN TOTAL COST(THEORETICAL) 


PROBLEM 

NO. 

7. IMPROVE- 
MENT OF 

HI 

7.IMPR0VE- 
ENT OF 

H2 

7. IMPROVE- 
MENT OF 

H3 

1 

11.88 

6.48 

16.76 

2 

21 .38 

19.71 

25 . 80 

3 

24.00 

67.43 

12 . 46 

4 

22.04 

21.29 

22 . 85 

5 

20.58 

12.11 

18.77 

6 

19.98 

-0.26 

19.98 

7 

32.33 

48.23 

65.04 

8 

14.84 

9 . 414 

20 . 79 

MEAN 

20 . 88 

23.05 

32.81 

S.D. 

6. 10 

23.11 

22.43 






A summary of results of the simulation carried out with eight dirferent 
problems are given in the tables 7.2.1 to 7.2.8. These tables contain the 
starting solution, the improved solutions of these heuristics with the 
percentage improvement. Simple mean and standard deviation of the 
percentage improvement of every heuristic on all performance measures 
considered was calculated and tabled in 7.2.9 to 7.2.12. 

7.3. Interpretation of Results 


The results obtained in this study gives a surprise. Particularly 

heuristic 1 improves the overall performance in a significant way. 
standard deviation of the improvement of heuristic 1 and 3 seem to be too 
much , but, the fact is that when some jobs get delayed and some jobs get 
expedited, it is very difficult to assess the dynamic impact on " shop 
dynamics ". 

On an average, the heuristic 1 improved the aggregate late ness by 
6.61%, actual inventory costs by 48.53 %, number of stock outs by 3.57 7 . and 
total cost by 20.68%. The heuristic 2 doesn’t improve the aggregate 
lateness in a significant way. But, on an average it improved actual 
inventory carrying cost by 34.39%, number of stock outs by 0.94 %, aggregate 
lateness by 0.88 % and total cost by 23.05%. The results are justified by 
the fact that when lot sizes are small, the aggregate number of components 


waiting f or other components will reduce which in turn reduce the actual 
inventory carrying costs. But, it may not be the same for aggregate 

lateness. Because, required amount of quantity may not be available for 

assembly. 

The heuristic 3 improves the actual inventory carrying cost in a 

significant way, because it does the things namely, decrease in the lot size 



and also changing the criticality. On an average this heuristic improved 
the actual inventory carrying cost by 54.60%. In the same time it increased 
the setup costs also. The improvement in the total cost using this 
heuristic is more than the increase in the total cost using heuristic 2. In 
aggregate lateness, the improvement is slightly less than that of heuristic 
1, which implies that there is no need of changing the lot sizes to reduce 
the aggregate lateness. In number of stock outs, the average improvement is 
3.58% which is equal to that of heuristic 1. 

Since, heuristic 3 over performs heuristic 2 in all aspects except the 
number of stock outs, heuristic 3 is prefered to heuristic 1. 



CHAPTER 8, CONCLUSION 


Previous research in Material Requirement Planning has been concerned 
largely with the single pass lot sizing and also with different dispatching 
rules in a job shop environment. Little work has been done to address the 
uncoordinated arrival pattern of different components, which is a common 
problem in the multi-stage, multi-product production system. As ailready 
stated, coordination is vital since shortage of a single item can halt the 
assembly process, resulting in sharply diminishing productivity and 
increased inventory costs. It is indeed a significant achievement, if it 
is possible to reduce the work-in-process inventory by even small but 
substantial percentage. 

With this in our mind, we developed a heuristic to schedule the 
arrival of sub-assemblies at an assembly centre in a more co-ordinated 
manner. Preliminary results indicate that such a heuristic improves the 
performance in a powerful manner. This heuristic could be refined in many 
ways to improve the performance of a dynaonic sjrstem such as MRP. 
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SOFTWARE DEVELOPED 



APPENDIX A 


SOFTHAKE DEVELOPED 



program mrpCinput, output); 


max_mc_no 

« 15; 

max_j)rocess_no 

= 5; 

max_assly_no 

= 45; 

max_child_assly 

* 5; 

max_period_no 

= 40; 

infinity 

= 32000; 

max_statisticsjpt_no 

= 4; 


type 


max_assly_mc_no = 16.. 20; 


event_type = (collect_statistics, 
processing gets over. 
assly_gets_over): 


ee s= record 
no, 

early_,demand 
waitin g t ime 
qty_needed 
arrival_time 
qty_arrived 
lot_number 
end; 


: integer; 

: real; 

:array [l..max_period_no) of integer; 
rarray ll..max_period_no] of read; 
:array [l..max_period_noI of integer; 
:array [l..max_period_no] of integer; 


ptr_to_lot = "lot; 


lot = record 

assly_no, 

size, 

current_mc_no : integer; 
operation_nuniber : integer; 
totad_wait_time, 
work_rem, 
due time. 



crr_value, 

Pt. 

tiine_of_joining_the_queue : real; 
next : ptr_to_lot; 
end; 

ptr_to_assly * ^ iot2; 

lot2 = record 

assly_no. 
current_mc_no, 
early_assly_order, 
order_completed, 
no_of_childs : integer; 

late_ness :real; 

child_infor :array[l..max_chiid_assly] of ee; 

next : ptr_to_assly; 

end; 

process_data *= record 

mc_no : integer; 

setup_time, 

processing;_time_j>er_unit, 
setup_cost : read 
end; 

child_assly = record 

child_assly_no, 

qty_to_be_contributed : integer; 
end; 

assly “ record 

assly_no, 

transf ormed_assly_no, 
parent_no, 
qty_for_parent, 
no_of_child_asslys. 



**o_of_machines_needed, 

•®c_no_for_assly, 

no_of_processing_needed : integer; 


cost,assembly_time: real; 


needsjprocessing : boolean: 

• * 

array_of _processing_sequence 

: arraytl..max_process_no] of process_data; 


array_of_child_asslys 

: arTay[l..max_child_assly] 
of child_assly; 


end; 

processing_machine = record 

work_injprogress : boolean; 

utilisation_time, 

finish_time : real; 

no_of _jobs_waiting tinteger; 

lot_ptr, 

Dtr to lot being processed : 
end; 


ptr_to_lot; 


assly_machine = record 

assly no being assembled. 
current_lot_size tinteger; 
work_in_progress : boolean; 
finish_time treal; 
assly_ptr :ptr_to_assly; 
end; 


(• newly declared types *) 


data_for_j)eriod * record 



demandl : integer; 
setup_cost, 

inv_carry_cost : real; 
end; 


tables *5 record 
scheduled_rec, 
projected_on_hand, 
gross_req, 
net_req, 
lot_quantities, 
planned_order : integer; 
ordering_cost. 
inv_carry_cost :real; 
end; 

inv_record *record 
qty_on_hand, 
open_reple_order, 
lead_time : integer; 
cost:real; 

aa :arrayIl..max_period_no] of tables; 
end; 

event_record = record 

event ; event_type; 
relevant_no : integer; 
time : real; 
end; 

var 

clock_time, 

simulation_time, 

waiting time during processing. 

initial_time, 

ordi^size, 

set_up__cost, 

inven_cost, 



total_cost, 

next_event_at, 

umi, 

end_late, 

best_late, 

cost ; real; 

continue, iterations_can_continue : boolean; 
temp_record :data_for_period; 
ooo,kk,mn,mmm,g,f,h,o,oo:text; 

input_data_for_ww :an'ayIl..niax_period_nol of data_for_period; 
iteration_required,best_iteration,no_of_stockout,no_of_setup, 
heuristic, stock_out,option,iteration_no, 
aaa,yy,current_period,i,no_of_period,mm,lot_no :integer; 
array_for_basic_data_about_asslys :array[L.max_assly_no] of assly; 
demand :array[l..max_assly_no,l..max_period_no] of integer; 
array_for_process_mc:array[l..max_mc_nol of processing machine; 
arcay_for_assly_mc :arTayImax_assly_mc_nol of assly_machine; 
item,base :ptr_to_lot; 
succ,pred :ptr_to_assly; 

inv_array : array! l..niax_assly_no] of inv_record; 
priority :arrayll..inax_assly_nol of integer; 
late_data :array[l..max_assly_no] of real ; 
aggre_late ;array[L.20] of real; 
inv_status :array[l..max_assly_nol of integer; 
k,past_history:array[l..max_assly_no] of integer; 
diff,there_cost,actual_cost :array[l..max_assly_no] of real ; 
r,s :set of O..max_assly_no; 

procedure read_assly_data: 

var 

ff :text; 
assembly_no, 
child_assly_no, 
needsjprocess, 
no,product_no, 
process_no 


:integer; 



begin (* read_assly_data •) 
assign(f f / asslyS. dat’ ); 

reset(ff); 
assembly_no :* 0 ; 
while not eof(ff) do 
begin 

assembly_no :=assembly_no + 1 ; 

writeln(assembly_no); 

readlniff, 

array_for_basic_data_about_asslys[assembly_no].assly_no, 

®^*'®JLf®r_basic_data_about_asslys[assembly_no]. 

no_of_child_asslys, 

needs_process, 

array_for_basic_data_about_asslys[assembly_no].cost); 
writelni array_f or_basic_data_about_asslys[ assenibly_no ] . assly_no ) ; 
array_for_basic_data_about_asslys[assembly_no]. 
needs_processing := (needs_process = 1); 

if (arTay_for_basic_data_about_asslys[assembly_no]. 
no_of_child_asslys > 0 ) 

then 

begin 

for i :=1 to array_for_basic_data_about_asslys[asseinbly_no]. 
no_of_child_asslys do 

with arTay_for_basic_data_about_£isslys[assembly_no] do 
with arTay_of_child_asslys[i] do 
readln(ff,child_assly_no,qty_to_be_contributed); 
end; 

if array_for_basic_data_about_asslys[assembly_no]. 

no_of_child_asslys > 1 then 
begin 

readln(ff,arTay_for_basic_data_about_asslys[assembly_no]. 

mc_no_for_assly. 

array_for_basic_data_about_assIys[assembiy_nol. 

assembly_time); 


end; 



If (array _f or_basic_data_about__asslysIassenibly_no3. 


then 


needs_processing) 


begin 


with ®rray_for_basic_data_about_assiys[assembly_nol do 
begin 


readin(ff,transf<MTned_assiy_no, 

no of processing needed): 

parent_no := 1; 
end; 


for process_no := 1 to 

arTay_for_basic_data_about_asslys[assembly_nol. 
no of processing needed do 

with array_for_basic_data_about_asslysIassembly_no] do 
with arTay_of__processing_sequence[process_no3 do 
begin 

readin(ff,inc_no,setup_time, 

processing_time_per_unit,setup_cost); 

end; 


end 


eise (* coiiect_parent_data •) 

with array_for_basic_data_about_asslys[assembly_no] do 
readln(ff,parent_no,qty_for_parent); 

end; 


close(ff); 

end; {* read_assly_data •) 
procedure read_assiy_demands; 

var 

input_demand_data :text; 
assly_no. 


aaa. 



period_no, 

no_of_demands, 

demand_no rinteger; 

begin (* read__assly_deTnands •) 

for assly_no :« l to inax_assly_no do 
for period_no := l to max_period_no do 
demandlassly_no,period_no] := 0 ; 

assigni input_demand_data, ’demands, dat ' ) ; 
reset(input_demand_data); 

while not eof (input_demand_data) do 
begin 

readln(input_demand_data,assly_no,no_of_demands); 

writeln(assly_no); 

for demand_no := 1 to no_of_demands do 
begin 

read(input_demand_data,period_no); 

readln{input_demand_data,demand[assly_no,period_no]); 

end; 

end; 

(• writelnC FROM DEMAND DATA FILE ’);•) 

for assly_no := 1 to max_assly_no do 
for period_no :* 1 to max_period_no do 
inv_array[assly_nol.aalperiod_no].gross_req := 

demandlassly_no,period_nol; 


close( input_demand_data) ; 
end; (* read_assly_demands •) 
procedure read_input_data; 


begin (• read_input_data *) 



read_assly__(lata: 

read_assly_demands; 

end; (* read_input_data •) 

procedure initialisation; 
var 

ii,nnn: integer; 
inv_data:text; 

begin 

assign(inv_data, ’inv8.dat’ ); 

reset(inv_data); 

for i := 1 to max_assly_no do 

with inv_array[il do 

begin 

writeln(i); 

(* readln(inv_data,qty_on_hand,open_reple_order,cost,lead_time); ♦) 

readln(inv_data,nnn,lead_time); 

writeln(nnn); 

op)en_reple_order := 0 ; 

qty_on_hand :* 0; 

for ii :*= 1 to max_period_no do 

with aa[ii] do 

begin 

read(inv__data,ordering_cost); 
orderin g cost := ordering_cost * 5 ; 

(• if array_for_basic_data_about_asslys[i].no_of_child_asslys > 1 then 

inv_carry_cost := 3.0 else inv_carry_cost := 1.0; *) 

scheduled_rec :*= 0; 

projected_on_h 2 “»d :* 0; 

net_req :*= 0; 

lot^quantities := 0; 

planned_order :* 0; 

end; 

for ii :*= 1 to max_period_no do 
with aalii] do 



begin 

read(inv_data,inv_carry cost); 
inv_carry_cost inv_carry_cost; 
end; 
end; 

for i:* 1 to max_assly_no do 

inv_arraytil.aall].scheduled_rec := inv_array[i].open_reple_order; 
close(inv_data); 
end; 

procedure do_lot_si2ing(n:integer); 

var 

optimum_policy :arTay[l..max_period_no] of integer; 

tenip_array :array[l..max_period_no,l..max_period_no] of real ; 

t_star,sum,ij,t,l,h,k,x,ii:integer; 

suml.minl,min2 :real: 

ff, order :array[l..max_period_no]of real; 

begin (* do_lot_sizing *) 

{* initialisation *} 

for i ;= 1 to max_period_no do 
with input_data_for_ww[il do 
begin 

demandl:= inv_arrayln).aali].net_req: 
setup_cost := inv_array[n].aa[i].orderingjcost; 
inv_carry_cost := inv_arTaytn].aa[i].inv_carTy_cost; 

end; 

(* for i :*= 1 to inax_period_no do writeln{input_data_for_ww[i].in'v_carry_cost 


t_star := 1; 
j :*= 1; 

1 :* 0 ; 


for i :* 1 to inaxjperiod_no do 



begin 

temp_arrayli,l] :• o.O; 
temp_arTayII,2] ;»0.0; 
end; 

for current_period :« l to maxj)eriod_no do 
begin 

with input_data_for_ww(current_period) do 

begin 
i :« 0; 

for j :»1 to curTent_period-l do 
begin 
suml :*0; 

if((current_period-j)*inv_carTy_cost • demandl < setup_cost ) then 
begin 

(• summation over h fc k •) 

for h j to current_period -1 do 
for k :* h+1 to curTent_period do 
suml :* suml +input_data_for_ww[h].inv_carry_cost 
• input_data_for_wwIk].demandl ; 


(• STORING IN A ARRAY •) 


1 := i+l; 

if j ol then temp_arrayll,l] := suml+iiq>ut_data_for_wwIjl.setup_cost 

+ ff[j-ll 


else temp_arTayIl,l] := suml 

+uq>ut_data_f<H'_ww[jl.iMti^_cost ; 


temp^arrayd.Zl := j ; 
end; 

end; (* for j< current period loop «ids *) 


{* calculation for the current period t *} 


1 :« 1 + 1 ; 

if current jperiod ■ 1 thMi teiiq>_anra3rll,l] seti^_cost 

else temp_arrayll,l] j® setiq)_cost + 
ff[current_period-ll; 



temp_arrayU, 2 ] ;■ current^period; 
end; (• WITH LOOP ENDS •) 

{• calculation of the minimum in temp.array •> 

mini temp_arrayll,i]; 

min 2 :■ temp_array[l, 2 ]; 

for i :« 1 to 1 do 

begin 

If temp_array[i,i] < mini then 
begin 

mini :*temp_arrayli,l]; 
min2 :* temp_arrayli,2]; 
end; 
end; 

ffIcurrent_periodl :« mini ; 
t_star :* round(min2) ; 
optimum_policy[curTent_periodl ;*= t_star; 
end; (• FOR ENDS •) 

{• finding the optimal Orders *} 

for i :■ 1 to max_period_no do orderli] := 0; 

i :« max _j)eriod_no; 

total_cost :* total_cost + ff[niax_period_no]; 
there_costIn] :* ffImax_period_nol; 

(• writein(ooo,n:5,’ theretical cost’,there_costtn]:8:2); *3 

mini :« there_cost[nl; 

repeat 

X :* optimum_policyIiJ; 
sum :*0; 

for ii :* round(x) to i do sum :® sum +input_data_for_wwliil.d«nandl; 
orderixl :»sum; 

mini :» mini - input_data_for_wwIxl.setup_cost; 
inv_array[nl.aa[x].lot_quantities := sum; 

i :»x-l; 
imtil X »1; 

lnven__cost :» inven_cost + mini ; 

set__up_cost :* sct_iq)_cost + ( there_costtn] “ mini ); 

(* writeln(ooo, 'setup cost *,there_costln] -minl:8:2); •) 



end; (• LOT SIZING ENDS *) 


procedure expIore(assly_no: integer); 
var 

j,cc: integer; 
begin 

with ^>'ray_for_basic_data_about_assIys[assIy_no] do 

if no_of__child_asslys > 0 then 

begin 

for j := 1 to no_of_child_asslys do 
begin 

with array_of_child_asslyslj] do 
begin 

cc:=array_for_basic_data_about_assIys[assly_no]. 
array_of_child_asslys[ j ]. child_assIy_no; 
for i := 1 to maxjperiod_no do 

inv_array[cc].aa[i].gross_req := inv_array[assly_no].aa[il.gross_req 
* arTay_of _child_asslys I j ] . qty_to_be_contr ibuted 

+ inv_atTay[cc].aa[il.gross_req; 

array_f or_basic_data_about_asslys[assly_no].parent_no := -1; 

array_f or_basic_data_about_asslysIcc].parent_no := 

array_for_basic_data_about_asslyslcc).parent_no - 1 ; 

end;{ with} 

end; 

end 

else array_for_basic_data_about_asslys[assly_no].parent_no := -1; 

end; 

procedure do_four_basic_steps(assIy_no: integer ); 

var 

j , cc, pp, order_time: integer ; 

begin 

(• writeln(f,assiy_no); 

writetf.’TimePeriod current’); 

for j := 1 to max_period_no do writeif,' ’,j:4); 



writeInCf); 

write(f, ’external dem ’); 

for j := 1 to max_period_no do write(f,’ ’.demand[assly_no,j]:4); 
writeln(f); 

writelf, ’gross •); 

for j := 1 to max_period_no do write(f,’ ’, 

inv_arraylassly_no]. aal j].gross_req: 4) ; 

writeln(f); •) 

(• netting •) 

with inv_array[assly_no] do 
begin 

(* write(f,’schd rec ’); 

for i := 1 to max_period_no do write(f,’ ’,aa[i].scheduled_rec;4); 
writeln(f); *) 

for i := 1 to max_period_no do 

with aa[i] do 

begin 

case i of 

1 : begin 

pp := qty_on_hand + scheduled_rec - 
inv_array[assly_nol.aa[i].gross_req ; 
if pp > 0 then 
begin 

projected_on_hand := pp; 
net_req := 0; 
end 
else 
begin 

projected_on_hand := abs(pp); 
net_req := abs(pp); 
end; 

end; 

2 ..inaxjperiod_no : begin 

if aali-ll.net_req = 0 then pp := 



aa[i-l].projected_on_hand + scheduled_rec 
inv_array I assly_no ] . aal i ] . gross_req 
else pp := scheduled__rec - 

inv_arrayIassly_no].aa[i].gross_req; 

if pp > 0 then 

begin 

projected_on_hand ;= pp; 
net_req := O; 
end 
else 
begin 

if aa[i-l].net_req = 0 then projected_on_hand := abs(pp) 

else projected_on_hand := 

aa[ i -1 ] . pro jected_on_hand +abs( pp ) ; 

net_req := abs(pp); 
end; 

end; 

end; { case} 
end; { for i > 

(* write(f,’proj hand ’,qty_on_hand:5,’ ’); 

for i := 1 to niax_period_no do writeCf,’ ’,aa[i].projected_on_h 2 uid: 4 ); 
writeln(f); 

write(f,’net req ’); 

for i := 1 to max_period_no do writelf,’ ’.aa[i].net_req:4); 
writeln(f); •) 
end; { with inv_rec } 

end; (• do _f our_basic_steps •) 

procedure mrp_initialise; 
var 

assly_no, ii, order_time .-integer; 

begin 

repeat 

for assly_no := 1 to max_assly_no do 

with array_for_basic_data_about_asslys[assly_no] do 



begin 

wrlteln(assly_„o;5,- -.parentjo); 

if ( par.irt_no = 0 ) then explore(aKl,_„o), 
end; 

ii := 1; 

for i := 1 to max_assly_no do 
begin 

if ( a’'ray_for_basic_data_about_asslys[i].parent_no <> 

then ii := 0 ; 
end 

until ( ii = 1 ); 
set_up_cost := 0.0; 
inven_cost := 0.0; 
total_cost := 0.0; 
no_of_setup := 0; 

for assly_no := 1 to max_assly_no do 
begin 

writeln(f); 

writeln(' ASSEMBLY NO ’,assly_no); 
do_f our_basic_steps(assIy_no); 
do_l ot_sizing( assly_no ) ; 

(* write(f,’lot ’); 

for i ;= 1 to maxjperiod_no do writeif,’ 
inv_array [assly_no]. aa[i]. lot_quantities: 4); 
writeln(f); •) 

for i := 1 to maxjperiod_no do 

if inv_array[assly_no].aa[i].lot_quantities > 0 then 

begin 

order_time := i - inv_array[assly_nol.lead_time; 
if order_time < 1 then order_time := 1 ; 
inv_array[ assly_no }. aa[order_t ime ] . planned_order : 

inv_array[assly_nol. aali]. lot_quantities + 

inv_array[assly_no].aa[order_time].pianned_order ; 

end; 

write(f, ’planned order 

for i := 1 to max_period_no do write(f,’ 

inv_array [ assly_no ] . aa[i I . pianned_order : 4 ) ; 



writeln(f); 


writeln(f,’Iead time 
end; 


, inv_array [ assly_no ] . lead_time: 4 ) ; 


for assly_no := 1 to max_assly_no do 
for i := 1 to max_period_no do 

if inv_arraylassly_no].aaIi].planned_order > 0 

then no_of_setup ;= no_of_setup + i ; 


end; { mrp_inio } 


procedure insert_in_assly_mc(item2:ptr_to_assly); 

var 

pred_ptr,succ__ptr :ptr_to_assly; 
begin 

with array_for_assly_mcIitem2".current_mc_no] do 

begin 

if assly_ptr = nil then begin 

assly_ptr := item2; 
end 

else 

begin 

predjptr := assly__ptr; 
succjptr := assly_ptr''.next; 
while succ_ptr <> nil do 
begin 

pred_ptr := succ_ptr; 
succ_ptr := succ_ptr^.next; 
end; 

pred_ptr'".next := itemZ; 
end; { else } 

pred__ptr := array_for_assly_mc[item2''.current_mc_no].assly_ptr; 
while(not{predjptr = nil)) do 
begin 

pred__ptr := pred_ptr^.next; 
end; { while } 
end; < with > 
end; { insert over } 



procedure load_the_mc(ij;integer;time:real); 

var 

pred_ptr, least, succ_ptr :ptr_to_lot; 

min_crr:real; 

pred_ptrl.succ_ptrl :ptr_to_assly; 
begin 

case ij of 
l..max_mc_no : 

with array_for_process_mc[ij] do 

if (( lot_ptr <> nil ) and ( work_in_progress = false ) 
and ( no_of_jobs_waiting > 0 ) 
and ( ptr_to_lot_being_processed = nil )) then 

begin 

mm:=mm+l; 

predjptr := lot_ptr; 
succ_ptr := nil; 
least := nil; 
min_crr := 100000.00; 

while pred_ptr <> nil do 
begin 

pred__ptr*.crr_value := 

(( pred_ptr''.due_time - time) - pred_ptr".work_rem ) / 
pred_ptr''.work_rem ; 

if predjptr^. crr_value < min_crr then 

begin 

array_f or_process_mc[ i j ]. ptr_to_lot_being_j)rocessed : = 
least := succ_ptr; 

niin_crr ••= pred_ptr^cr^_value; 

end; 

succ_ptr := pred_ptr; 


pred_ptr; 



predjptr := pred_ptr",next; 


end; 

if { least a nil ) then 

lot_ptr := lot_ptr^.next 

else least^.next := ptr_ to lot being processed'^.next: 

array_for_process_mc[ij].ptr_to_lot being processed^. next := nil; 

work_in_progress:= true; 
no_of_jobs_waiting := no_of _jobs_waiting - 1; 

with array_for_basic__data_about_asslys[ptr to lot bein g p rocessed''. 

assly_nol do 

with array, of processing sequencelptr to_lot_bein g p rocessed^. 
operation_number] do 

array_for_process_mc[ij].finish_tinie ;= time + 
ptr to lot being processed'.pt: 

if finish_time > simulation_time then 

utilisation_time := utilisation_time + 

( simulation_time -time ) 

else utilisation_time := utilisation_time + 

ptr to lot being processed", pt; 

with ptr to lot being processed^ do 
total_wait_time := total_wait_time + 

( time - time_of_joining_the_queue ); 

base := array_for_process_mc[ij].lot_ptr; 

while base <> nil do 
begin 

base := base^.next; 

end; 

end; 
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end; { case } 

end; { for all assl } 

procedure insert(itemme:ptr_to_lot); 

var 

pred_ptrl.succ_ptrl :ptr_to_lot; 
ii: integer; 

begin 

with &>'ray_for_process_mcIitenune''.current_mc_no] do 

begin 

base ;= lot_ptr; 
if lotjptr = nil then begin 

lot_ptr := itemme; 
end 

else 

begin 

pred_ptrl := lot_ptr; 
succjptrl := lot_ptr'‘.next; f 
while succ__ptrl <> nil do '* 

begin 

pred__ptrl := succ_ptrl; 
succ_ptrl := succ_ptrr.next; 
end; 

pred_ptrr.next := itemme; 
end; { else } 

no_of _jobs_waiting:=no_of_jobs_waiting + 1; 

pred_ptrl := array_for_process_mcIitemme".current_mc_no].lot_ptr; 
while(not(pred_ptrl = nil)) do 
begin 

predjptrl := pred_ptrr.next; 
end; { while } 


end; { with) 



end; {insert } 


procedure load_the_assly_mc(timee:real); 

var 

aa,iij,m rinteger; 
w_t :real; 

load,ini_found ; boolean; 
predjptr :ptr_to_assly; 

wait_time :array [l..max_child_assly] of real ; 
begin 

for ii := 1 to max_child_assly do 
wait_time[ii] ;= 0.0; 

for ii := 16 to 20 do 

with array_for_assly_mclii] do 

begin 

(* writelnCass me no ’,ii); •) 
predjptr ;= assly_ptr; 
j:=0; 

while predjDtr <> nil do 
begin 

with pred_ptr" do 
begin 

load := true ; 

for j := 1 to no_of_childs do 
begin 

if early assly_order >= child__inforf j].early_demand 
then load := false; 

end; 

case load of 
true : begin 



If cun-cmj>=rl.d > «rl,.assl,.„rd.r then 
no_er_stockout no.of_stoek(,ut t 1; 

with arra,_for_assl,.mc|predj>u-.en,ient.mn no] do 
begin 

currentJot_si2e := 

in v_array I pred_ptr . assly_no ] . aa[ early_assly_order ] . 
planned_order; 

work_in_progress := true ; 
finish_time ;= timee + 

array_for_basic_data_about_asslysIpred_ptr".assly_no]. 

assembly_tiine; 

assly_no_being_assembled := pred_ptr*.assly_no: 

if finish_time > early_assly_order * 4000.00 then 
late_ness := late_ness + 

(( finish_time - early_assly_order * 4000.00 ) * 
current_lot_size ) ; 

order_completed := order_completed + 1 ; 

for j := 1 to no_of_childs do 
with child_infor[j] do 

begin 
m := 1 ; 
aa := 0; 

with array_f or_basic_data_about_asslys[pred_ptr''. 

assly_no] do 

with array_of_child_asslysIj] do 

if child_assly_no = child_infor[j].no then 

begin 

aa := array_of_child_asslys[j].qty_to_be_contributed 
in v_ar ray [pr ed_ptr assly_no ] . aa[ ear ly_assly_order ] . 
planned_order; 
end; 

while ( chiIdJnfor[j].qty_arrivedlin] < 1 ) 

and ( m < 40 ) do m := m + 1; 



repeat 


if childJnfor[j].qty_arrived[m] > o then 

begin 

if child Jnfor[j].qty_arrived[inl <= aa then 
begin 

wait_time[j] qty_arrived[m] 

( timee - arTival_time[m] ) + wait_time[j]; 
aa := aa - qty_arrivedlm]; 
qty_arrived[m] := -1 ; 
end 

else 

begin 

wait_time[j] := aa 

* ( timee - arrival_time[m] ) + wait_time[j]; 
qty_arrived[m] := qty_arrived[m] - aa ; 
aa := 0; 
end; 

end; 

m:= m +1; 

until ( aa = 0 ) or ( m = 41 ); 

wait_time[j] := wait_time[j] / 4000.00; 
end; { for j } 

for j := 1 to no_of_childs do 
begin 

aa 0; 
w_t := 0.0; 

with array_f or_basic_data_about_asslysIpred_ptr''.assly_n 

with array_of_child_asslys[jl do 
if child_assly_no = child_infor[j].no then 
begin 

aa := array_of_child_asslys[j].qty_to_be_contributed • 
inv_array[pred_ptr".assly_no].aa[early_assly_order].plan 



end; 


inv_status[childjnfor[j].no] := 
inv_status(child_infortj].no] - aa ; 

w_t := aa • wait_time[j] • 

inv_arraylchildJnforlj|.„ol.aa[cmTent_i«rIodl.i„v_c»r 

:= W_t / currentJot_size ; 

child_infor(j].waiting_time := child_infor[j]. waiting ti 

end; 


child_infor[j].qty_neededIearly_assly_order] := -i; 

j:=early_assly_order + 1 ; 
ini_found := false ; 
repeat 

if ( inv_arraylassly_no].aa[j].planned_order > 0 ) then 
begin 

early_assly_order := j ; 
ini_found := true ; 
end; 

j := j +1: 

until (( ini_found = true ) or ( j = max_period_no + 1)); 

if ( j = max_period_no + 1 ) then 

begin 

early_assly_order := 1; 
end; 


end; 

end; 


end; < case > 


end; < predjJtr } 



pred_ptr :» pred_ptr''.next; 


end; { pred <> nil } 
end; { with array for all me } 
end; { load the assly me } 

procedure ®®^'^^'_fo>'_parent(tt,amt:integer;times:real); 
var 

np.nq,nn,ii.j,jj,m: integer; 
predjptr.succjptr :ptr_to_assly; 
over, load.f ound, ini_f ound : boolean; 

begin 

lot_no := lot_no + 1 ; 

case array__for_basic_data_about_asslys[ttl.parent_no of 

0 : begin 

inv_status(tt] :*= inv_status[tt] - amt ; 
end; 

1 : begin 

nn 0; 
nq :== 0; 

ii :*= 16; 
found :» false; 
over :* false ; 

while { ii <* 20 ) and ( found = false ) do 
begin 

with array„for_assly_mc[ii] do 

begin 

predjptr :* assly jptr; 

j;*0; 

whil= ( prsdjtr o ril ) and ( found = falao ) do 



begin 


with predjDtr" do 
begin 

for j :*= 1 to no_of_childs do 
begin 

if child_infor[j].no = tt then 
begin 

found := true ; 

jj := child_infor[j].early_demand; 
chi ld_inf or[ j ] . qty_arrived [ current_period ] 


repeat 

jj :« child_inforIj],early_demand: 

if ( amt >= childJnforIj].qty_needed[jj] ) and 
( childJnfor(j].qty_neededIjj] > 0 ) then 

begin 

amt := amt - child_infor[j].qty_needed[jj] ; 
child_inforIj].qty_neededIjj] := 0; 

writelnCh, ’fully satisfied’); 

m jj +1 ; 
ini_found := false ; 
repeat 

if ( child_infor[j].qty_needed[m] > 0 ) then 
begin 

child_infor[j].early_demand := m ; 
ini found :® true ; 


end; 



m := m +1; 


until (( ini_found 


true ) or ( m = max_period_no + 1 )); 


end 


else 

begin 

childjnf orl j].qty_needed[ j j] ;= 

childJnfor[j].qty_needed[Jj] - amt ; 
amt := 0; 
end; 


until ( amt 0 ) or ( m = meCK_period_no ) or ( over = true ); 

childjnf or(jl.arrival_time[current_period] := times; 

childjnf or[j].lot_numberlcurrent_period] := iot_no ; 
end; { if child found } 

end; { for all childs } 

end; { with pred } 

predjptr :» predjitr^.next; 

end; < while pred <> nil } 

end; { with me ii } 

ii := ii + i ; 

end; < while ii < 21 } 

end; { case parent > 
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begin 


while ( amt > 0 ) do 
begin 

nn :* maxjperiod_no; 
nq := 0; 


for ii := 16 to 20 do 

with array_for_assly_mc[iil do 
begin 

(• writelnCass me no 
assly_ptr; 

J:*0; 

while pred_ptr <> nil do 
begin 

with pred__ptr" do 
begin 

(• writeln(’assly_no*,assly_no);*) 
for j :» 1 to no_of_childs do 
begin 

(• writelnCchild ’,child_infor[j].no:3);*) 
if child_infor[j].no = tt then 

begin 

( • writelnlh, ’child’, child_inf or[ j].no, 

’early’ , child_inf or[ j]. early_demand, 

’qty needed ’,child_infor[j].qty_needed[child_inforIj]. 
early^demand]);*) 

if child_infor[jl.early_demand < nn then 
begin 

nn :« child_inforIj].early_demand; 
nq :« child_inforIj).qty_needed[nn]: 
end 



else 


if ( childJnforIj].early_demand = nn ) and 

child_inf orl j]. qty_needed[child_inf or[ jl.early_demand] 

< nq ) then 

nq ^i'*i*i_infor[j].qty_neededInn] ; 
end; { if child found ) 

end; < all child } 

end; { with prd ptr } 

pred_j)tr :* pred_ptr''.next; 

end; { while pred <> nil } 

end; { for all me > 

{ early demand for all parents found } 

{ going inside and satisfying early demand } 

writeln(h, ’after searh nn ’,nn:4,’nq ’,nq:4); 

if nq <> 0 then 

for ii :» 16 to 20 do 

with array_for_assly_mc[iil do 

begin 

(“ writelnCass me no ’,U): *) 
pred__ptr :* asslyjptr; 
j:»=0; 

while pred _ptr <> nil do 
begin 

with predjptr' do 


begin 



for j :« 1 to no_of_chilcls do 
begin 

(• writelnCchild ’,childJnfor[j].no:3);«) 
if ( ( child_infor[j].no = tt ) and 

( chilciJnfor[jl.early_demand = nn ) and 
( child_infor[j].qty_needed[nn] = nq ) ) then 


begin 

jj := child_infor[j].early_demand; 

(• writelnih/ parent no :’,assly_no,’ early period 

child^inforl j].early_deinand,’ qty needed 

child_infor( j].qty_needed[child_infor[j].ean'Iy_demand]:5); *) 


if ( amt < nq ) then 
begin 

child_infor(J].arrival_time[current_period] := times; 

(• writeIn(’c’,child_infor[j].qty_needed[current_period]); *) 
child_infor[jl.qty_arrived[current_period] := amt + 
chi Id^inf orl j ) . qty_arrived( current jperiod ] ; 
child_infor(j].qty_needed(nn] := child__infor[j].qty_needed[nn] 

“ amt ; 

child_infor[jl.lot_number[current_period] ;= lot_no; 

(• writeln(’cc’,child_inforIjl.qty_arrived[current_period]); •) 

amt :« 0; 
end 
else 
begin 

child_infor[j].arrival_timeIcurrent_periodl := times; 

childJnfor[j].qty_arrived[current_periodl := nq + 

child Jnfor(jl.qty_arrived[current_periodl ; 

child_infor(J].lot_number[current_periodl := lot_no ; 
child JnforIjl.qty_needed[nn] :=0; 

amt :» amt - nq ; 

(• writeln{h,’cc’,child_infor[j]. 

qty_arrivedlcurrent_period]); *) 



m :* nn +1 ; 
ini__found := false ; 
repeat 

if ( child JnforIj].qty_needed[ml > 0 ) then 

begin 

childJnfor[j].early_demand := m ; 

ini_found := true ; 
end; 

m := m +1; 

until (( ini_found = true ) or ( m = niax_period_no + 1 )); 

end; 

{• writelnih, ’assly no ’,assly_no,’anit remaining ’,amt); *) 
end; (• if child is found •) 

end; 

{ for all child } 

end; < with pred ptr > 
pred_ptr :« predjptr".next; 

end; { while pred ptr <> nil } 

end { with me } 

else begin 

(• writeln(oo,’some problem ’); *) 
amt :« 0 ; { take care of this } 

end: 

■{ most demanding person satisfied } 

{♦ writeln(h/amt after %amt:4); •) 


end; < while amt > 0 } 



end { if parent > 2 } 


end; { case ends } 


load_thc„assly_mc( times); 


end; { search for parent ends } 

procedure initialise_the_newjtem{i,ss:integer;time:real); 

var 

j , pr oces;;„no, i i ,yy : integer; 

frt)nt:ptr_to_lot; 

item?,: ptr_to_assly; 

begin 

(• writelniin’ ASSEMBLY NO initiated with ’,i, ’size ’,ss:6) 
inv_statusnj :* inv_status[i| + ss ; 

case array_for_basic_data_about_asslys[i].needs_processing of 

true : 

begin 

new(item): 
witli item" do 
begin 

assly_no :« i; 

size :® ss; (* inv planned order was there *) 

toial_wait_time 0.0; 

crr„value :*= 0.0; 

if heuristic * 1 then 

due_time :« current jperiod * ordi_size + 

( past„history[iJ • 0.1 * ordCsize ) 
else due^time :* currentjeriod * ordi_size ; 

time^of Joining„the_queue :« 0.0; 
operation^number :*=1: 



with «rra,_for.batic_data_ab«ut_a.sl,sll) do ' 

begin 

With arTay_„r_procaaato^a.,da„ca|op.,«,„„ 

begin 

current^mc^no :* me no; 

pl:«setup_time + processing_time_per_unit * size; 

end; 

work_rem:«0.0; 

for process_no := i to no_of_processing_needed do 
with array_of_processing_sequence[process_no] do 

work_rem :« work_ren,+ setup_time processing_time_per_unit • ss 
{• writelnCh, ’work_remaining’,work_rem:6:2);*) 

{•crr_vaiue :* 4000.00 /work_rem; 
writeln(h,'CRR VALUE ’,crr_value:4:2):*) 
end; {• array •) 


time_of Joining;_the_queue := time; 
next :« nii; 
end;(* item •) 

insert (item); 

if array^f orjprocess_mc[item''.current_mc_nol. 
work_injprogress * false then 

load_the_mc( item^. current_mc_no,clock_time) ; 


end; 

false : search^for jparentii.ss.time); 

end; (• case for needs processing end *) 

end; (• for ail asslys *) 

procedure update„the_item(item:ptr_to_lot;time:real); 

var 

predj>ir,succjptr sptr^to^assly; 
j.tt,ii,ss: integer; 
mnj :real; 



begin 
il :’=2; 
with item'' do 
l>egin 

work_rem :« work^rem -pt; 

itcm'.operation^nuinber := item^operation_number +1; 
witii array^for^basic_data_about_asslys [assly_no] do 
if item^operation_numbe^ <= array_for_basic_data_about_asslys [assly_ 
>f'o_oi'jP>'ocessing_needed then 

begin 

witii nrr ay__of _jprocessing_sequence[operation_number] do 

t>egin 

current_mc_no :* mc_no; 

pt :» setup_,tiine + processing_time_per_unit • size; 
end; 
ii:«l; 

end {then ) 

else 

begin 

inv_,status[item''.assly_no] := 
inv_statuslitem''.assly_no] - item'.size ; 
tt :»= array_for_basic_data_about_asslys[assly_no]. 
transf ormed„assly_no; 

il:* 0; 
end; { else } 

end ;{ Item) 

case ii of 
0 ; begin 

ss :* 
j :» I 


item'*.size; 



while . 

J :» J +1 : 

mnj :« time - ( j + i) • 4000 .O; 

if mnj > 0.0 then late_data[ttl := late_data[tt] + 
mnj • item'', size; 
waitlng^.limc_duringjprocessing 
i t cm'' . iotal_wait_time: 


ss :* item", size; 
dispose(item); 

initialise_the„new_item(tt,ss,time); 

end; 

1 : imgin 

itcm''.tlme„orjQining_the_queue :*= time ; 
insert (item); 

if array„forjprocess_mclitem".current_mc_no].work_in_progress = fal 
ioad^t he_mc( item". current^mc_no,clock_time) ; 

end; 

2 

end; 

end; { update ti>e item ) 


procedure sequencing initialisatiQns(init_time:real); 


var 

j, au, yy . process^no, ii , m ; integer; 

front:plr_to_iot; 

item2, pred jptr : ptr_to_assly , 

begin (• »equencingL_initiaUsations *) 

(• wr ilelnOi/SEOUENCING INIT AT ’,imOime:5;2); •) 


for i I to IS do 



while ^^v_ar^ay(item^assIy_nol.aa[j].pla^ned order . 
J J +1 ; 

rnnj time - { j + i) • 4000.O; 

if mnj > 0.0 then late_datalttl := late_data[tt] + 
mnj • item". size; 
wiUting^time^duringjproGcssing :« 
item''.total_wait_time; 


ss :« item", size; 
dispose! item); 

i n i t ia 1 i se_t he_ne w Jtem(tt , ss, time ); 

end; 

! : begin 

i!em".time„ofJoining_the_queue := time ; 
insert (item); 

if «»rr-iiy„forjproccss_mc(item".current_mc_no].work_in_progress = fal 
load_the„mc(item".current_mc_no,clock_time); 

end; 

2 

end; 

end; { iiptlate the item ) 


procedure seutiencing initiaUsations(init_time:real); 

var 

J, aa, y y . pr0(;ej;s„no, i i , m : Integer; 

front :ptr.„to_ lot; 
i t em2 , pred j;>t r ; pt r_to_assly ; 

begin (• sequencing^initialisations •) 

(• writclidh/SKOUKNCINC INIT AT ’,imtjime:5:2): *) 


for i 1 to 15 do 



while 

j : = 


inv^arraylltem", 
j +1 : 


assly_no).aaIj].pianned_order = 


mnj :« time - ( j + i) . 4000.0; 
if mnj > 0.0 then late_dataltt]’ :» 

mnj • ^lem^size; 

i>"e-.durlng_j)rocessing ;« 
itenr.toial_wait_time: 


late_data[tt] + 


ss do 


;« item''. size; 
dispose(item); 

initiaUse_the_newJtem(tt,ss,time); 

end; 

1 : begin 

itemMime^ofJoininohe^queue := time ; 
insrrtiitem); 

if *»•I•^‘y^for•„process_mc(item^cu^^ent_mc_no].wo^k_in_p^og^ess 
load , t hr^ mc( it enr. current_mc_no,clock_time); 

end; 


end; < update the item } 


pitjcetlure sequencing„initialisations(init_time:real); 

var’ 

j , aa . y y , pr( »rcss. no, i i , m : integer; 

frnnt:pti 

itemi*,pre<l,.ptr:ptr_tc)_assly; 
begin (• jiequencing_initialisations •) 

(• wrilelnCh/SKQUENCING INIT AT ’,init_time:5:2); *) 


fal 


fw i 1 t« 15 du 



with array_for_process_mcIi] do 
if finish_time > init_time then 
begin 

if finish_time < simulation_time then 
utilisation_time := finish_time - init_time 

else utilisation_time := simulation_time - init_time ; 

end 

else utilisation_time := 0.0; 

for i := 1 to 15 do 

with array_for_process_mc[i] do 

if ( finish_time < init_tiine ) and { finish_time <> 0.0 ) then 
begin 

finish_time := 0.0; 
ptr to lot being processed := nil ; 
work_in_progress := false; 
end; 

yy := 0 ; 

for i := 1 to max_assIy_no do 

begin 

case au~ray_for_basic_data_about_asslys[i].no_of_child_asslys of 
0 : 

if inv_array[i].aa[current_period].planned_order > 0 then 

begin 

initialise_the_new_item(i, inv_arrayl i ]. aa[ current_period] . 
planned_order, init_time ) ; 
yy := yy + 1; 
end; 

1:; (* writeln(h,’it has only one child ’);*) 

2..max_child_assly :; 


end; (* case for no of child *) 
end; 

end;(* sequencing over *) 



procedure start_thejprocess; 
begin 

for i := 1 to max_mc_no do begin 

load_the_nic(i,0.0): 

end; 

end;{ stau't_the process } 


procedure event(tiine:real); 
var 

order,n,m rinteger; 
init_over :booleein: 

begin 

ni:=l; 

for i:= 1 to max_mc_no do 

with array_f or_process_mc[i] do 

if Dtr to lot being processed <> nil then 

begin 

if time = finish_time then 
begin 
m := 0 ; 

item := ptr to lot being processed: 
Dtr to lot being processed := nil; 
work_in_progress := false ; 
finish_time := 0.0; 
load_the_mc(i,time); 
update_the_item( item, time) ; 
end; 
end; 

if m = 1 then 

for i := 16 to 20 do 

with array_for_assly_mc[i] do 

if finish_time = time then 

begin 



end; 

procedure 

var 


(* writeln(oo,assly no being assembled. 

assly oprn over’);*) 
finish_time := 0.0; 
work_in_progress := false; 

init_over := false; 
order := 0; 
n := 1; 
repeat 

if demand! assly no being assembled, n] > 0 then 
begin 

order := n; 
init_over := true ; 
end; 

n:= n+1; 

until ( init_over = true ) or ( n = current_period + 
if ( ( order > 0 ) and ( init_over = true )) then 
begin 

current_lot_size := current_lot_size - 

demandfasslv no being assembled, order ] ; 

demand! asslv no being assembled.order] := 0; 
end; 

initialise the new itemCasslv no being assembled, 
curr ent_l ot_size, time ) ; 

end; 


update_ciock_time(time: real ) ; 


init_over : boolean; 
nn:real; 



begin {• update_clock_time •) 
i:=l; 
nn:= 0; 

next_event_at := 0.0 ; 
init_over:»=f alse; 
repeat 

if array_for_process_mc[i].finish_time > 0.0 then 
begin 

nn := array_for_jprocess_mc[i].finish_time; 

init_over := true; 

end; 

i:= i+1; 

until (init_over) or (i > max_mc_no); 

if init_over = f alse then 
repeat 

if array_for_assly_mc[i].finish_time > 0.0 then 
begin 

nn := array_for_assly_nicIi].finish_time; 

init_over := true; 

end; 

i:= i+1; 

until (init_over) or (i > 20); 

if ( nn = 0.0 ) then clock_time := simulation_time+100 
else 
begin 

(• writelnC’init event at :’,nn:5:2); *) 
next_event_at := nn; 

for i := 1 to max_inc_no do 
begin 

if ( array_forjprocess_mcli].finish_time > time ) and 
(array_for_process_mc[i].finish_time < next_event_at) then 
next_event_at := array_for_process_mcIi].finish_time; 
end; 



for i := 16 to 20 do 
begin 

if (array_for_assly_mc[iI.finish_time > time) and 
{array_for_assly_mcti].finish_time < next_event_at) then 
next_event_at := array_for_assly_mc[i].finish_time; 
end; 

clock_time := next_event_at; 

end; 

end; (• update_clock_time *) 


procedure do_sequencing; 
var 

check_time :real; 
itemZ :ptr_to_assly; 
j,m :integer; 
ini_found :boolean; 

begin (* do_sequencing *) 

(* initiadisation of array_for_process_mc *) 

for i := 1 to max_mc_no do 
with array_for_process_mc[i] do 
begin 

work_in_progress := false; 
utilisation_time := 0.0; 

finish_time := 0.0; 
no_of_jobs_waiting := 0; 
lotjptr := nil; 

ptr to lot being processed := nil; 
end; 

for i := 16 to 20 do 

with array_for_assly_mc[i] do 



begin 

current_lot_size := 0; 
asslv no being assembled 0; 
asslyjptr :» nil; 
finish_time :»= 0.0; 
work_in_progress :»=false; 
end; 

lot_no := 0; 

current_jperiod := 1; 

for i := 1 to max_assly_no do 

if array_for_basic_data_about_asslys[i].no_of_child_assIys > 1 then 

begin 

new(item2): 

with array_for_basic_data_about_asslys[i] do 

with itemZ'' do 

begin 

iate_ness := 0.0; 

order_conipleted := 0; 
current_mc_no := mc_no_f or_assly; 

assiy_no := array_for_basic_data_about_asslysli].assly_no; 
no_of_childs := no_of_child_asslys; 

j:=l; 

ini_found := false ; 
repeat 

if ( inv_array[assly_no].aa[j].planned_order > 0 ) then 
begin 

early_assly_order := j ; 
ini_found := true ; 
end; 

j := j +1; 


until {( ini_found = true ) or ( j = maxjperiod_no + 1 ) ); 



for j :w 1 to no_of_childs do 
begin 

child_inforlj).no :* array_of_child_asslysI jl. 

child__assly_no; 

child inf or[ il. waiting time := 0.0; 

for m := 1 to max_j)eriod_no do 
begin 

child_infor[jl.arrival_time(in] := 0.0; 

child_infor[j).qty_arrived[m] 0; 

child_infor[j].qty_neededIml :* 0; 

child^inforlj].lot_numberfml :* 0; 

if inv_array[assly_no].aaIm].planned„order > 0 then 
child_infor(j].qty__needed(m] := 
array_of _child_asslysl j ] . qty_to__b®_co>nit>'ibuted 
* inv__array [assly_no ) .aaf m ] . planned_order 

else child_inforIjl.qty^needed(ml -1 ; 


end; 


m :»= 1; 

ini_found :* false ; 
repeat 

if ( child_inforlj].qty_neededlm] > 0 ) then 
begin 

child_inforlj].earIy_demand := m ; 
ini_found := true ; 
end; 

m := m +1; 

until ( ini_found *= true ) or ( m > max_period_no ); 


end; 


next nil; 



end; { item > 


insert_ln_assly_mc(item2); 

end; < if it is a assembly } 

clock_time := 0.0; 
continue := false; 


while ( current_period < 41 ) do 
begin 

(• writelnC' CURRENT PERIOD :’,current_period); *) 
initial_time := (current_period - 1 ) * 4000.00; 
simulation_time := current_period * 4000.00; 

sequencing initialisationsC initial_time ) ; 

for i ;= 1 to max_mc_no do load_the_mc(i,initial_time); 
update_clock_time( initial_time) ; 
while(clock_time <= simulation_time) do 
begin 

event(clock_time); 
update_clock_time(clock_time) ; 
end; 

(• writelnimn, ’current period’, currentjperiod:5); 

writelnimn,’ me no utiIisation_rate ’,’jobs w 

for i := 1 to max_mc_no do 
with array_for_process_mc[il do 

writelnimn,’ ’.i:4.’ ’.utilisation_time/4000. 

* ’,finish_time:8:2,’ ’,no_of_jobs_waiting:5 



for i := 1 to max_assly_no do 
begin 

(• writeln(inv_status[i]); *) 
cost := cost + ( inv_status[i] * 

in v_array [ i 1 . aa[ current_period ] . inv_carry_cost ) ; 

actual_cost[i] := actual_costlil + ( inv_status[i] * 

inv_array [ i 1 . aa[ current_period 1 . inv_carry_cost ) ; 

end; 

(* writeln(’cost’,cost:10:2); *) 

current_period := current_period + 1; 
continue := false ; 

for i := 1 to max_inc_no do 

if array_for_process_nic[i].utilisation_time < 4000.00 then 
umi := umi + array_for_process_mc[i].utilisation_tinie/4000.00 
end; { while } 

end; (* do_sequencing ends *) 

procedure collect_and_print_statistics; 

var 

predjptr :ptr_to_assly; 
j.yy.yyy.cc anteger; 
min, max :real; 

begin (* collect_andjprint_statistics *) 
writelni 00 , ’ statistics’ ) ; 

for i := 1 to max_assly_no do priorityli] := 0; 

for i := 16 to 20 do 

with array_for_assly_mc[i] do 



begin 

(• writelnCass me no •) 

pred_ptr :»= assly_ptr; 

J:*0; 

while predjptr <> nil do 
begin 

with pred_ptr" do 
begin 

(* writelniooo.’ass no ’,pred_ptr".assly_no); *) 
min := 1000000000.00; 
max := 0.0 ; 
yy := 0 ; 

yyy == O; 

for j := 1 to no_of_childs do 
begin 

if child inf or[ il. waitin g t ime > meoc then 
begin 

yy := child_infor[ jl.no ; 
max := child inf or[ il. waiting time ; 
end; 

if child inf or [ il. waiting time < min then 
begin 

yyy := child_infor[ jl.no ; 
min := child inf or [il. waiting time ; 
end; 

end; < all child over } 

stock_out := stock_out + order_completed ; 

case pred_ptr'‘.assly_no of 
1,5,19.27 : begin 

if order_completed <> 0 then 



begin 

late_ness /order_coinpleted/l 00000: 8:2); 
end_late :* end_late + ( late_ness /order_completed); 
end else end_late := end_late + late_ness ; 
end; 

2. .4,6. . 18,20. .26,28. .max_assly_no : ; 

end; 

late_datalpred__ptr^.assly_nol := late_ness ; 

if yy <> 0 then 
begin 

prioritylyy] := prioritylyy] + 1 ; 

with array_for_basic_data_about_assIys[yy] do 

if no_of_child_asslys > 0 then 

begin 

for j := 1 to no_of__child_asslys do 
begin 

with array_of_child_asslyslj] do 
begin 

cc: *:array_f or_basic_data_about_asslys[ assly_no ]. 
array_of _child_asslys[ j ]. child_assly_no; 
prioritylcc] := prioritylcc] + 1 ; 

end; 
end; 
end; 
end; 

if yyy <> 0 then 
begin 

prioritylyyy] := priority[yyyl - 1; 

with array_for_basic_data_about_asslyslyyy] do 

if no_of_child_asslys > 0 then 

begin 

for j := 1 to no_of_child_asslys do 
begin 

with array_of_child_asslys[j) do 
begin 



cc:»«array_for_basic_data_about_asslys[assly_no]. 
array_of _chi ld__assly s[ J ]. chi ld_assly_no; 
priority! cc] :* priority [cc] - 1 ; 

end; 
end; 
end; 
end; 

end; 

pred_ptr ;= predjptr''.next; 
end; 
end; 

for i := 1 to max_assly_no do 

with array_for_basic_data_about_asslys[i) do 

if no_of_child_asslys > 1 then 

aggre_lateliteration_no] :*= aggre_late(iteration_no) + late_data[j] 

umi := umi / ( max__period_no * max^nic_no ); 
for i := 1 to max_assly_no do 
case priority! i] of 
0 :; 

1..50:past_history!il := past__history!i] + 1 ; 

-50.. -1 : past_history!i] := past_history!i) -1; 

end; { case } 

if ( heuristic = 2 ) and ( iteration_no > 1 ) then 
if best_iate > aggre_lateIiteration_no] then 
begin 

best__late := aggre_late!iteration_no] ; 

best_iteration := iteration_no; 

end; 

end; (• collect_and_j)rint_statistics *) 
procedure do_heuristic_two; 


var 



: integer; 


i , iii , old_lots, no, new_lots 
smt,ratio,max :real; 
stop : boolean- 

old :arrayIl..maxjperiod_no] of integer; 

begin 

r :» [ ]; 

for i := 1 to max_assly_no do 
begin 

diffii] := ( actual_cost[i] ) - there_cost[l] ; 

if (not(i in s ) ) then 

begin 

ratio difflil / rostli! ; 

if ( ratio >0.1 ) then r + I 

end; 

end; 

max :* 0.00; 
no :* 0; 

for i :* 1 to max_assly_no do 

if ( i in r ) and (not(i in s )) then 

if diff[i] > max then begin 

max :* diffli); 
no i ; 

end; 

stop := false ; 
iii := 0; 
repeat 

iii := iii + 1; 

for i := 1 to maxjperiod_no do inv_array[nol.aa[i].inv_carry_cost 
inv_arraylno].aali].inv_carry_cost • 1.?. ; 

old_Iots := 0 ; 
new_lots:*=0; 
smt :* 0; 

for i := 1 to maxjperiod_no do 
begin 

if inv_arraylnol.aaIi].lot_quantities > 0 then 



begin 

old_lots := old_lots + 1; 

smt := smt + inv arrav[no].aa[i]. ordering cost: 
end; 

if inv_array[no].aa[i].net_req > 0 then new_lots := new_lots + 1; 

old[i] := inv_array[no].aa[i].lot_quamtities ; 

end; 

no_of_setup := no_of_setup - old_lots ; 
inven_cost := inven_cost - ( there_cost[no] - smt ); 
set_up_cost := set_up_cost - smt ; 
total_cost := total_cost - there_cost[no]; 
writelnC ’old’ ,old_lots, ’total net’ ,new_lots); 

old_lots := new_lots: 
do_lot_sizing(no ) ; 
new_lots := 0; 

for i := 1 to max_period_no do 

if inv_array[no].aati].lot_quantities > 0 then new_lots := new_lots + 1; 
for i := 1 to max__period_no do 

if oldli] <> inv_array[no].aa[i].lot_quantities then stop := true ; 
writelni ’new’ ,new_lots); 
no_of_setup := no_of_setup + new_lots; 
until ( stop = true ) or ( new_lots >= old_lots ); 
k[nol := k[no] + 1 ; 
if k[no] >= 2 then s := s + [no]; 


procedure mrp_system; 
begin 

ordi_size := 4000.00; 
iteration_no := 0; 

for i := 1 to max_assly_no do 
begin 

priority] i] := 0; 

past_history[i] := 0; 

late_data[i] := 0.0; 

actual cost[i] ;= 0.0; 



klil := 0; 
end; 

for i :*= 1 to 20 do aggre_late[i] := 0.0 ; 

r := I 1: 

s := [ ]; 

repeat 

umi := 0.0 ; 

cost :=0.0; 

stock_out := 0; 

end_late := 0.0 ; 

no_of_stockout: =0; 

for i := 1 to max_assly_no do 

begin 

late_data[i] := 0.0; 

actual_cost[i] := 0.0; 

diff[il := 0.0; 
end; 

iteration_no := iteration_no + 1; 
writeln( ’ITERATION NO :’,iteration_no); 
writeln(ooo, ’ITERATION NO :’,iteration_no); 

iterations_can_continue :=fzJse; 
while not(iterations_can_continue) do 
begin 

do_sequencing; 

iterations_can_continue := true; 
end; 

collect_and_print_statistics; 

if heuristic = 2 then do_heuristic_two; 

{ disposing all asslys } 

for i := 16 to 20 do 

with array_for_assly_mc[il do 



if assly_ptr <> nil then 
begin 

pred :>= asslyjptr; 
repeat 

succ := pred”'. next; 
if pred <> nil then dispose(pred); 
pred := succ; 
until ( succ = nil ); 
end; 

for i := 1 to max_mc_no do 
with array_for_process_mc[i] do 
if lot_ptr <> nil then 
begin 

item := lot_ptr; 
repeat 

base := item", next; 
dispose(item); 
item := base; 
until ( base = nil ); 
end; 

until ( iteration_no = iteration_required ); 
if heuristic = 2 then writeln(ooo,’best sol late : 

’iteration no 


begin (* main *) 

assign{f,’ol6’); 

rewrite(f); 

assign(h,’ss’); 

rewrite(h); 

assignioo, ’result.dat’ ); 
rewrite(oo); 

assign! mn , ’ resultZ. dat ’ ) ; 
rewrite(mn); 


’,best_late:10:2, 

’ ,best_iteration: 4) ; 



assign(kk, 'result3.dat’ ); 
rewrite(kk); 

assign(ooo, ’result4.dat’ ); 
rewrite(ooo); 
read_inpt‘ data; 
initialisation; 

write( ’ENTER CHIOCE OF HEURISTIC NO (1/2/3) :’) 
read(heuristic); 

mrp_initialise; 

read_assly_data; 

best_late := 100000000000.00; 
best_iteration := 0; 
case heuristic of 
1,2 : begin 

iteration_required := 20; 
mrp_system; 
end; 

3 : begin 

heuristic := 2 ; 

iteration_required := 20; 

mrp_system; 

initialisation; 

read_input_data; 

mrp_initialise; 

read_assly_data; 

iteration_required := best_iteration; 
heuristic := 2; 
inrp_system; 
heuristic := 1; 
iteration_required := 20; 
inrp_system; 
end; 

end; 

close(f); 

close(h); 

closeCoo); 



close(mn); 

close(kk); 

close(ooo); 



